DIFFUSIVITY OF LATTICE GASES 



JEREMY QUASTEL AND BENEDEK VALKO 

Abstract. We consider one component lattice gases with a local dynamics and a stationary prod- 
uct Bernoulli measure. We give upper and lower bounds on the diffusivity at an equilibrium point 
depending on the dimension and the local behavior of the macroscopic flux function. We show 
that if the model is expected to be diffusive, it is indeed diffusive, and, if it is expected to be 
superdiffusive, it is indeed superdiffusive. 



1. Introduction 

In this article we consider lattice gas, or particle models for the stochastically forced diffusion 
equation with non-linear drift in M. d 

d t p = V -j(p) + V -DVp+VDV (1) 

where j(p) is the macroscopic current, D is a diffusivity and £j, i = 1, ... ,d are smooth processes 
approximating independent space-time white noises. Our main interest is on the long time, large 
wavelength behaviour of solutions, in particular of the space-time correlation functions in equilib- 
rium, and how this depends on the non-linearity j(p) and the dimension d. 

In one dimension only, ([!]) coincides with the stochastic Hamilton-Jacobi equation for a height 
function h defined through V/i = p, 

d t h = j{Vh) + DAh + y/D£. (2) 

The case j(x) = v\x\ 2 with white noise forcing is the Kardar-Parisi-Zhang (KPZ) equation. In one 
dimension, it has a deep connection with representation theory and integrable systems, which has 
recently led to some distributions being computed exactly [7] . Note however that the scaling of the 
stochastic Hamilton-Jacobi equations in higher dimensions are different from ([I]). Most importantly, 
with white noise forcing, have formal white noise invariant measures in all dimensions, while 
([2]) have them only in d = 1. 

To explain the heuristics for the scalings, we write an approximate equation for the rescaled 
variable 

Pe (t,x) = e- a p(e- fS t,e- 1 x) 
with a small parameter e representing the ratio of small to large scales, and a, f3 > 0, 

d tPe = e a ~P +1 e- a j' t (e a Pe)d XlPe + e^ +2 V • DVp € + e'^'^ VDV ■ £. (3) 
i=l 

The approximation is in the last term which we scaled as if it had correlations E[£(t, x)£(s, y)] = 
5{t — s)5(x — y). In this case of Gaussian white noise, the rescaled distribution (,(at, bx) would be 
statistically the same as a , X . 

1 
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Now one has several natural choices for a and /3 depending on the dimension d, the nondinearity, 
as well as what it is that one desires to see. 

Hyperbolic (Euler) scaling (/3 = 1). If we take a = we obtain 

d tPe = V • j(p e ) + eV • DV Pt + VDV ■ £, (4) 

from which it is not hard to guess that p t —> p satisfying 

d tP = V-j(p). (5) 

Weak solutions of ([5J are not unique. However they are, for bounded po, [361 US] if supplemented 
by the entropy condition, for all c £ M, in the weak sense, 

dt\p - c| < V • sgn(p - c)(j(p) - j(c)). (6) 

It is expected that the Euler scaling of asymmetric interacting particle models leads to the entropy 
solutions in wide generality, but it has only been verified in special cases [261 EE1 Ell CEZ1 CO 

Longer time scales (/3 > 1). If we look on longer time scales, the scaling now depends on the 
behaviour of j(p) near po = 0, as well as on the dimension d. We will expand j around po = 0, and 
note that in the resulting Taylor series, the first two terms j(0) and j'(0)p can be easily removed 
from the equation by a simple change of variables and a coordinate shift. So we set j(0) = j'(0) = 0. 
The first non-trivial term in the series is therefore ^j"(0)p 2 . Generically, one has j"(0) ^ 0, in 
which case j'(e a p t ) ~ e a p e j"(0). This means that we can take (3 = 2 (the diffusive scale) and a = 1 
to get, 

d 

8 tPt = e- a j'i(e a Pe)d XiPe + V • DVp t + VDV ■ £. (7) 

i=l 

In d = 1, the noise dominates, and clearly there is no limit, d = 2 is critical, and one expects 
logarithmic divergences. Only in d > 3 is there a limit |§1 121]. 

d t p = V-j"(0)p + V-DVp. (8) 

So we can obtain a limit by diffusive scaling only if d > 3. 

To see the fluctuations, one is guided by the natural rescaling of the invariant white noise, which 
leads one to guess that a should be d/2. We arrive at 

dtp e = e^-P it ^ d/2 j'^ d/2 Pe )d XiPe + 6 2 ^V • DVp e + e^VDV- £. (9) 
i=l 

In d > 3, the clear choice is = 2, and the fluctuations are given by an (infinite dimensional) 
Ornstein-Uhlenbeck process [El H7] 

d tP = V ■ DVp + v^DV • f . (10) 

In d = 1, we are constrained by the first term. In the generic case j"(0) ^ 0, we are forced to 
take f3 = 3/2. This is the KPZ scaling. Note however that the limit process is not the naive guess, 
which would be the entropy solution of 

d t p = \j"{v)d xP 2 (n) 

because it can be checked by the Lax-Oleinik formula that this does not preserve (standard) white 
noise [10] . while ^ does, for each e > 0. In fact, the limit process is a Burgers equation forced by 
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a highly non-trivial, and not well understood residual noise. Some finite dimensional distributions 
are known, but not all [8]. 

One can see that d = 2 is critical for Q , and as usual in the critical case one expects logarithmic 
corrections. On the other hand, if we are at an inflection point, then the behaviour is different. If 
if (0) = 0, then e~ a j' i (e a p e ) is of order e a if j"'(0) ^ and of order e 2a if j"'(0) = 0. One then checks 
that the above d > 3 scalings now apply in d = 2 if all j"(0) = and in d = 1 if j"(0) = j"'(0) = 0. 
These are the diffusive cases. The final case is that of a generic inflection point in d = 1. We have 
j"(0) = but j'"(0) / 0. Referring to (Jok, we have e~ d / 2 f (e d / 2 p e ) ~ e d l 2 and we can see once again 
that (3 = 2 but we are in a critical case and therefore there should be logarithmic corrections. 

The model cases of ([I]) are thus j(p) = v ■ p n with n = 2 the generic stochastic Burgers equation, 
diffusive in d > 3, logarithmically super-diffusive in d = 2 and KPZ super-diffusive in d = 1; n = 3 
is the generic inflection case, diffusive in <i > 2, and logarithmically super-diffusive in d = 1; and 
ra > 4 is the double inflection case, diffusive in all dimensions. 

But these arguments are at best heuristic, because ^ is at best approximate, relying on the 
exact rescaling of pure white noise forcing £. However, no existence or uniqueness theorem in known 
for ([!]) with white noise forcing for non-trivial j(p) except for the very special d = 1, j(p) = p 2 
case, and even here one is restricted to a finite interval with periodic boundary conditions |14j . 
In fact, the expectation is that ([I]) does not have non-trivial solutions except in that special case. 
There are several ways around this. One could try to deal with the smoothed noise, however such a 
study is made difficult by the fact that we do not know the resulting invariant measures, which are 
used extensively in the analysis. Lattice gas models provide natural discretizations where one does 
know the invariant measures, so they are a natural place to start. As we will see, the superdiffusive 
behaviour is reflected in a diffusion coefficient which we can study. Our method is based on the 
Green-Kubo formula, which requires one to solve a resolvent equation in the space generated by 
the invariant measures, or, equivalently, to study space-time correlations of the flux. This requires 
detailed knowledge of the invariant measures. 

The superdiffusive scalings in the lattice gas models were identified at the physical level using 
the formalism of mode-mode coupling theory. Let p(k, t) be the spatial Fourier transform of the 
density p(t,x). S(k,t) = (p(—k,0)p(k,t)) is called the intermediate scattering function. The 
brackets denote expectation in equilibrium. An equation can be written for dfS(k,t), but, of 
course, it involves higher order correlations, equations for the time evolution of which involve yet 
higher order correlations, etc. The mode-mode coupling formalism approximates the higher order 
correlations at the first step by multilinear forms in S, producing a closed equation. A description 
of the procedure can be found in [15]. The result for ([!]) with j(p) = v ■ p n is 

d t S(k, t) = -D\k\ 2 S{k, t) - c(v ■ k) 2 {S(k, 0))- 1 [ S(k, t - s) S * ■ ■ -*S(k, s)ds (12) 



o 



n times 



where S indicates that this is only an approximate equation for S. The * denotes convolution in 
k, c is a constant. Assuming S(k,t) = S(0, Q^ e - D \ k \ 2t - a ( v - k ) 2t (^°& t )' : an d then solving for a and £ by 



identifying both sides of (12) as k — > oo one guesses that £ = 2/3 in the critical case n = 2, d = 2 
and C = 1/2 in the critical case n = 3, d = 1. These predictions were then backed up by extensive 
numerical simulation |34|. 
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Our results can be summarized the following way. For a large class of lattice gas models, which 
have Bernoulli invariant measures, if the model is expected to be diffusive, it is indeed diffusive, 
in a sense to be made precise in the next section. And, if it is expected to be superdiffusive, it is 
indeed superdiffusive. 

Recently, there has been a large amount of activity around such models in connection to the 
Kardar-Parisi-Zhang equation. For example, jT3j study the speed-changed exclusion models for the 
special d = 1, j(p) = cp 2 case of or Q. Following [5], they study the weakly asymmetric limit, 
proving that under that scaling the density fields are tight, and satisfy a plausible weak version 
of the Kardar-Parisi-Zhang equation, though, unfortunately, such solutions are not known to be 
unique. In comparison, our work focuses on the case of general flux j and dimension d, and we do 
not introduce small parameters into the problem, but study the scale of fluctuations directly. In 
this sense, our results are more in the spirit of, but far less exact, than the scaling results for the 
special solvable one dimesional models such as TASEP, PNG and q-TASEP. On the other hand, 
and this has to be emphasized, we are in a situation where there appears to be no exact solvabilty, 
and, except in the special d = 1, j(p) = cp 2 case, there is not even a known solvable model in the 
universality class. One can of course ask if some of these methods can be extended to models with 
more than one particle per site, such as zero-range, or where more than one particle moves at the 
same time. We have not pursued this here, but as long as one has the explicit invariant measures, 
the method should work, though it could be extremely complicated in practice. One can also study 
related continuous spin models |27j . For Gaussian invariant measures, they turn out to be very 
similar to the exclusion cases; the case of general potentials has not been studied yet. 

A finer question than the diffusivity would be the limiting fluctuation fields. As usual, what 
one expects is Gaussian fluctuations in the diffusive and logarithmically super-diffusive cases. The 
diffusive case is studied in [17]. The logarithmically super-diffusive case is still open. Only in 
the case of two dimensional exclusion with asymmetry in one direction are the scaling results fine 
enough that one might hope to prove the Gaussian fluctuation limit. In other cases, one still has 
to get the scale correct. 

2. Model and main results 

We now describe the model and results precisely. We consider lattice gas or speed changed ex- 
clusion models on Z d with local interaction and exclusion rule. They consist of particles performing 
continuous time random walks, with jump rate depending on the local configuration, and with the 
exclusion rule that jumps to occupied sites are suppressed. We think of it as a Markov process 
with state space is {0, l} z<i , the or 1 indicating the absence or presence of a particle at x G 
The infinitesimal generator is given by 

L f(v) = E E Kl/, T-zV) %(1 - Vx+y)(f(v X ' X+V ) ~ M) (13) 

x y 

where t x denotes the shift {r x r])y = r\ y - x -, and r : Z d x {0, l} zd — > M + gives the rate r(y,T- x r]) of a 
particle at x to jump to x + y. rj x (l — rj x+y ) indicates that there is a particle at x and no particle 
at x + y, so that the jump can be performed, and f(rj x ' x+y ) — f(rj) measures the change in the 
function / from the pre-jump configuration rj, to the post-jump configuration r] x ' x+y , which has the 
occupation variables at x and x + y exchanged. 
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We need to make some assumptions on the jump rates: 

Local evolution. For every y ^ the function r(y, •) > is a local function, i.e. only depends 
on the values rj z , \z\ < K where K is a fixed constant. Note that we can assume without loss of 
generality that r(y, rj) does not depend on (770, %)• Moreover, r(y, •) = if \y\ > K. 

Divergence conditional There exist local functions R±, . . . ,Rd so that 



Coercivity. The additive group generated by those y 6 Z rf for which r(y, rj) + r(—y, T y rj) > for 



Let us comment on the conditions. The locality of the rates is natural. It just means that 
a particle looks in a finite neighbourhood to determine its jump rate. The fact that under this 
condition there is a unique Markov process r/(£) on {0, 1} Z with e tL f(rj) = E[f(r](t)) | r/(0) = 77], 
follows from general theory [22 . 

Assuming we have local rates, the divergence condition implies that the product Bernoulli(/9) 
measures, ir p , p £ (0, 1) are invariant. This is our key assumption, which allows us to start 
the analysis. It is not a generic property, but on the other hand one can find a rich family of 
examples. The proof that it implies the invariance is as follows: By general theory the local 
functions form a core for the generator [22] , so to check the invariance one need only prove that for 
such functions / Lfdm p = 0. Using the fact that the change of variables 77 — > rj x ' y preserves tt p so 



The summation in x gives a telescoping sum and by choosing a large enough box the surviving 
terms from the sum lie well outside the support of /, and are therefore independent of it. Hence 
the last term vanishes. In fact, for local rates the other direction is also true; the invariance of 
the Bernoulli measures implies the divergence condition. This will not be used in our results; the 
proof of the statement is left to the interested reader. The divergence condition is the analogue in 
our context of the fact that the flow generated by a divergence free vector field preserves Lebesgue 
measure. 

The coercivity condition is not a strong assumption. It is equivalent to the ergodicity of the 
process and also equivalent to the symmetric part of the generator L being comparable to the 
generator 



d 



r (2/> */) ~ Vo) = ^e.Riiv) where V e J(t/) = f(rj) - f(r ei r]). 



(14) 




all 77 is equal to 7L d . 





(15) 



of the symmetric simple exclusion process. More precisely, for the symmetric part S 
exist < c\ < C2 < 00 such that the Dirichlet forms satisfy 



L+L* 
2 



there 



ci(f(S )f) < (f(-S)f) < c 2 (f(-So)f). 



(16) 



Sometimes referred to as the gradient condition [TS]. However, it is not the same as the gradient condition used 
in the theory of hydrodynamic scaling limits. 
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The Dirichlet forms are given by 

(/(-So)/) = !K p [£(/(^ +ei ) - f(v)) 2 } (17) 

and (f(S)f) = i^ p E^(r(y,r_ x .r ? ) + r(-y,r_ :E+?/ 7 ? )(/(r ? ^)-/(r ? )) 2 ]. We can perform rf*+» 
by doing ||y||i = m nearest neighbour switches. By the Cauchy-Schwarz inequality, [f{rj x,x+v ^ — 



f{v)) 2 — m Y^i=o (J '(jh' " 2j+1 ) — f(Vj)) 2 - Since r(-, •) is bounded, we get the right inequality of (16) 



The left inequality works similarly, but here we have to replace the switch iHilej using jumps 
from the set {y : r(y,rj) + r(—y,T y rj) > for all rf\. The coercivity condition makes this possible. 

We now give some explicit examples of rates. 

(1) Finite range exclusion process. Here r(y,rj) = p(y), this clearly satisfies the first two 
conditions. If we assume p(-) is such that the associated random walk can get to any site 
in Z d then the coercivity condition is satisfied as well. The special case p(l) = lu=i is the 
symmetric simple exclusion process. 

(2) The following rates provide a simple model satisfying the conditions in one dimension 
without being a finite range exclusion process 

K 1 ) V) = 3 - V-l - m, r(-l,r]) = 2, r (y,rj)=0 if \y\ ^ 1. (18) 

At certain points of the proof we will illustrate the computations using this simple model 
before proceeding with the general proof. 

(3) Let d = 1 and < y G Z. Then it is easy to check that both r{y, rj) = rji ■ ■ ■ f]y—\ and 
(1 — 771) ■••(!. — %-i) satisfy the divergence condition. (Similar products can also be defined 
in the y < case.) This gives the following model: suppose that we have an exclusion 
process satisfying the conditions, but we increase the jump rates of size y with a certain 
constant c y if there are no particles (or no holes) between the starting point and the end 
point. If {y : c y > 0} is finite then this model will satisfy all the conditions^ 

(4) It is easy to construct higher dimensional models from the one dimensional examples. Sup- 
pose we have one dimensional rate functions r\, ■ ■ ■ , which satisfy the conditions. Let 
ei, . . . , ed denote the usual unit coordinate vectors in Z rf and let P{r\ be the projection of 
i] G {0, 1} Z to its i th coordinate. Now define the (i-dimensional rate function r the following 
way: 

r(yei,r]) =ri(y,PiT]), y G Z and r(z,rj) = Q liz^ya. (19) 
It is not hard to check that this will satisfy our conditions. 



(5) The model ( 18 ) can also be extended into higher dimensions easily. Consider a d-dimensional 
exclusion with jump law p y , y G Z rf satisfying the conditions. Now modify the jump rates 
so that for a finite subset of y's the jump rate is changed to 

r(y,y) = p y {3 - rj-y ~ V2y) (20) 
The resulting model satisfies all our conditions. 



2 This is similar to, but not the same as the model q-TASEP, where r{y,rf) = unless y = 1, in which case 
r (!> rj) = Ys^oi 1 ~ 1 n )Vn+2 nri? 1 i 1 ~ r)i). The q-TASEP does not satisfy the diver gence condition (as well as having 
infinite range), and the invariant measures are not Bernoulli (see [6]). 
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We now define the flux. To define the microscopic flux we first write 

Lr lo = ( r (y> T yV)v-y(l - Vo) ~ r (v> v)vo(l - r] y )) (21) 
y 

Note that this is a finite sum and it can be written as a divergence J2i=i ^eiWi with local functions 
Wj. The vector W = (W\, . . . , Wd) is called the microscopic flux^ The flux for the adjoint L* is 
defined analogously; we denote it by W* = (Wf, . . . , Wjg). The macroscopic flux is 

j(p)=E p (Wx,...,W d ). (22) 

The entries of this vector are polynomials in p, 

3i(p) = Y2,y' e i E p r (y^ 7 l) r lo( l - %)• (23) 
y 



The model is asymmetric in the i th direction if ji(p) is not a constant. Note the abuse of language 
in that there are non-symmetric models with zero drift which are therefore not asymmetric. So we 
should really say "non-zero drift". However, the term asymmetric has become ubiquitous in the 
field and we will stick with it. We emphasize that in order to identify the flux, one needs to know 
at least something about the invariant measures. In fact, we do not know any model where the 
flux has been identified without knowing the invariant measures exactly. 

Examples. For a Id exclusion process with jump rate p(-) the microscopic flux is the linear 
combinations of terms of the form rj x (l — r] x + y ) and the macroscopic flux is bp{l — p) where b is the 
first moment of p(-). For the simple Id model defined in ( |18[ ) the microscopic flux is 170 ( 1 — ) ( 1 — 
— 77-2) and the macroscopic flux is j(p) = (1 — 2p)p{l — p). Note that this has an inflection 
point at p = 1/2. 

In general, one is interested in the large scale behaviour of the space-time correlations, 

S(x,t) = X ' 1 (ri(x, t); »/(0,0)> (24) 

where \{p) = /°(1 ~~ p) an d (/; g) is just a notation for the equilibrium covariance E[fg] — E[f]E[g]. 
Note that we will also use the notation (/, g) for the scalar product Efg. The expectations are 
with respect to ir p with some fixed p £ (0, 1). One easily checks [23] 

J2S(x,t) = l, and 52xS(x,t)=j'(p)t. (25) 

X X 

So the first non-trivial momenlQis the time dependent diffusivity 

D(t) = t-\- x (]T* 2 Mi), r,o(0) - pfj - (j'(p)t) 2 . (26) 

Note the normalization so that standard diffusive behaviour corresponds to D(t) = 0(1) and 
superdiffusive behaviour to D{t) — > 00 as t — > 00. For technical reasons our key observable will 



Note that although W is not unique, we can just fix a function which is a linear combination of terms of the form 

Note however, that only in the attractive case (e.g. exclusion) will one necessarily have S(x, t) > 0. In that case 
D(t)t is the variance of the position of the second class particle at time t, so it is actually connected to the second 
moment of a random variable. 
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be the Laplace transformed diffusivity given by the Green-Kubo formula (see Lemma 18 in the 
Appendix) , 



Du(X) := A 2 / V(x i -jKp)t) 2 5(x,i)e- A ^ = C ll + 2 X (K,(A-L)- 1 ^))-2 X ((^,(A-L)- 1 ^)), 
J o x 

(27) 

where C& = ^ y yf E r(y,vi), for local functions f,g, 



e-i i-\ 



■= E(/; r ^) = I™ 7 <£r i /;5><7> (28) 



1=0 «=0 



and, for Pi, the orthogonal projection to the linear subspace generated by the functions rji we define 
w i = ^{W i + W*-r 1 W i -ViW*), Vi = ^{Wi-W*). (29) 

Note that in the case of finite range exclusion processes, the third term on the right in equations 
(27) vanishes. This is because in that case V{ = V ei g with q being a first degree polynomial and 
thus the corresponding scalar product is zero. For the general case we will show that there is a 
model dependent constant C < oo such that 

(A — L)~ x vifj < C. (30) 

Hence our results are equivalent to statements about the blowup as A \ of ((wi, (A — L)~ 1 Wif). 

We are interested in the following predictions for the asymptotics for the diffusivity as A \ 0. 
We say that a model is diffusive (in the i th coordinate) if -Dm (A) has a finite limit and superdiffusive 
(in the i th coordinate) if Du(\) — > oo. The predictions make precise in this context the arguments 
of the introduction: 
Dimension d = 1 

(a) If j"(p) 7^ then the expected behavior is superdiffusive with D(\) ~ A -1 / 3 . 

(b) If j"(p) = 0,j"'(p) ^ then the expected behavior is superdiffusive with D(X) ~ (logA -1 ) 1 / 2 . 

(c) If j"(p) = j"'(p) = then the model is expected to be diffusive. 
Dimension d = 2 

(d) If j'{(p) then the expected behavior is superdiffusive with Dii(X) ~ (log A -1 ) 2 / 3 . 

(e) If j" (p) = then the model is expected to be diffusive. 

Dimension d > 3 The model is expected to be diffusive with D(X) converging to a finite limit. 

Our results say basically that if the model is expected to be diffusive, it is indeed diffusive, and 
if it is expected to be superdiffusive, it is indeed superdiffusive. Note that the constants in the 
statements may depend on the density p G (0, 1). 

Theorem 1. Assume that the speed change model has local, coercive rates satisfying the divergence 
condition. Then 
Dimension d = 1 

(a) lfj"(p)^0then 

dX~ 1/4 < D(X) < C2X- 1 / 2 (31) 

(b) If j"(p) = and f'(p) ^ then 

Ci log log X- 1 <D(X)< C 2 log A" 1 (32) 
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(c) Ifj"(p) = f"(p) = then £> U (A) < C. 
Dimension d = 2 

(a) If j'{(p) 7^ then we have 

Ci(log A" 1 ) 1 / 2 < £> U (A) < CalogA- 1 (33) 

W Iff{(p) = 0then D U (X)<C. 
Dimension d > 3 

In a// cases Da{\) < C . 

Remark 2. The d = 1, = 0, j"'(p) ^ case bears many similarities in scaling of the diffusivity 
to the d = 2 lattice gas models for Navier-Stokes studied in [20J and certain random diffusions in 
random environment in d = 2 |33j . This is not completely coincidental as there are some structural 
similarities between models in d = 1 at inflection points and models in d = 2 without preferred 
direction which make the both superdiffusivity and the estimates based on degree 3 test functions 
coincide in the two cases. However, the technical issues involved in the proofs are very different, 
and the test functions themselves are different, although the basic strategy based on variational 
methods is the same. 



The key issues in the present proofs are (i) the identification of the leading order term (see (42)), 
and (ii) to show that all the many other pieces of the asymmetric part of the generator do not play 
a role. In the generic d = 1, j"(p) 7^ case, for example, the leading order terms are the same as 



for asymmetric exclusion, and the eventual test function used for the lower bound (31) is the same 
as in earlier work [19] . What is new here in that case is (ii), i.e. the universality. 

In earlier work [241 125] , we proved the universality of the i 1 / 3 law among one dimensional exclusion 
processes with drift. This was possible because in that case we have a strong version of (ii), as 
well as a special solvable model (TASEP) within the class. Unfortunately, we are not aware of any 
solvable model in the d = 1, j"(p) = 0, j"'(p) ^ class, so at the present time we have no way 
to obtain the exact asymptotics (logi) 1 / 2 . Our lower bound is via a degree 3 test function, which 
is complicated, and in fact we only have identified its Fourier transform. In the d = 2, generic 
case of j"{p) 7^ for i = 1 or 2, one does have a type of exact result for the the special nearest 
neighbour model which is symmetric in one coordinate direction and totally asymmetric in the 
other coordinate direction, at p = 1/2 [37] . In principle, one expects that with sufficient effort, 
one could prove a strong enough result of type (ii) above, to show that this is universal among such 
models. We have not pursued this here. But note that even extending the results of [37] to p ^ 1/2 
is non-trivial as the diagonal part of the asymmetric part of the generator no longer vanishes, and 
is not straightforward to control. 

Note that the full predictions for the superdiffusivities can be obtained formally from the vari- 
ational method if one assumes certain scaling properties of the optimal test function, and (more 
severely) that all off diagonal terms in the computations vanish. The formal computations are 
described in [20] [271 [33]. Since they are fairly analogous in our case, we do not repeat them here. 

We note that there are many ways to measure diffusivity; we concentrate on D(\) here because 



our method is basically to study the blowup of the second term on the right hand side of (27). 
One can certainly study the blowup of the resolvent for functions other than the current, but 
the current seems most natural from our point of view, and through the Green-Kubo formula is 
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connected directly to the second moment of the correlation functions. From the Laplace transform 
to this second moment seems to us a technical issue about regularity, and it is hard to say whether 
one should prefer one to the other. In one dimension, one also has a height function h(t,x), 
whose discrete derivative is given by rj(t,x). Part (b) in Theorem [l] is roughly equivalent to the 
conjecture that if j"(p) = 0,j"'(p) ^ then h(t,Q) ~ ct + t 1 / 4 (logt) 1 / 8 C- The fluctuation ( is 
predicted to be Gaussian. Note that this is not related to anisotropic KPZ in d = 2, which is 
proved to have h(t,0) ~ ct + (logt) 1/,2 £ or isotropic KPZ in d = 2 which is conjectured to have 
h(t, 0) ~ ct + t a24 - '£. In order to study a fluctuation like t 1/ ' 4 (logt) 1 / 8 one would need a different 
method than the present one. The problem is that only objects whose time correlations decay 
slower than i -1 will have a diverging resolvent. Of course, a natural thing to try would be some 
quadratic form like (w, (A — L)~ 2 w) or perhaps a higher power (w, (A — L)~ p w), depending on the 
context. However, recall Lowner's theorem, which says that a function / preserves the definiteness 
of A — L if and only if it is real analytic and has analytic continuations to the upper and lower half 
planes which maps the upper half plane to itself. Among powers this singles out (w, (A — L)~ l w). 

To understand the microscopic origin of the conditions on j"(p) and j"'(p) we study the repre- 
sentation of functions in L2({0, l} d ,ir p ) in the following orthonormal basis. For any finite A C Z rf , 
let 

va=v p a = ru~ i/2 fe-p), ( 34 ) 

xeA 

where x = p(l — P)- Let M. n denote the local functions of degree n, i.e. / is a finite sum / = 

:iai= 



5jA|=n IaVA- We also denote by M. n the closure of this space in L 2 (P), i.e. 



feM n O / = E /A<7A, E fl < °°' ( 35 ) 

|A|=n |A|=n 



For f,g £ M n we have 



(f,g) = E Asa- (36) 

[A|=n 



If / G M n , 9 £ A4 m , m ^ n then we have (/, g) = 0, so our Hilbert space is naturally graded 

L 2 ({0,l} zd ,vr p ) = e~ A4 n . (37) 
For the ((, )) scalar product of /, g 6 M. n (defined in p8j )) we have 

((f,g))= E A5A+y = E Asa, /a = E A+y, (38) 

|A|=n,yeZ d AeB„[d] 



where B n [^] is the set of equivalence classes of size n subsets of Z rf modulo the translation. (I.e. 
A ~ A' if A' = A + y.) We call these dimension reduced variables. If d = 1 then we can represent an 
element of B n [l] as an element of Z" _1 : the equivalence class of the ordered ra-tuple (xi, x%, . . . , x n ) 
is represented by the vector (x2 — x\ — 1,2:3 — x-i — 1, . . . , x n — x„_i — 1). In that case / can be 
represented as a function on Z™ _1 and the scalar product is the usual: 

«/,<?» = (/,S>= E /»s(z)- (39) 
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Now suppose that (dropping the subscript i for clarity), 

™ = E<*?A (40) 
A 

is the representation of i th coordinate of the symmetrized microscopic flux. Note that ji(p) '■= E^pWi 
and ji(p) = E np Wi only differ in a linear term so their second and higher order derivatives will be 
the same. In fact, with a slight abuse of notation, we will drop the tilde in j from this point. Since 



E 



P+eVx 



X 



-1/2 



£ we get 



UP + e) = E <E p+ sV P a = E < £W x 



p_|A|_.-|A|/2 



and hence 



9 k P m = k\ x - k/2 e 



w 



A' 



|A|=fc 



Note that Wi is always a function of degree 2 and higher. Hence, 



and 



f{p) = o 



j"'(p) = o 



<3> 



E" 

|A|=2 

E'< 

|A|=3 







0. 



(41) 
(42) 

(43) 
(44) 



The problem becomes one of describing the behaviour of ((w, (A — L)^ 1 !!])) as A \ according to 
conditions such as the right hand sides of (43) and (44). 

We study ((w, (X—L)~ 1 w)} using the variational formula, which we learned from S.R.S. Varadhan, 

((w, (A - L)-M> = sup {2((w, /)) - ((/, (A - S)f}) - ((Af, (A - S^Af))} (45) 

where the supremum is taken over local functions. The operator A = ^ is the asymmetric part 
of the generator. An easy upper bound can be obtained by dropping the final, non-negative, term 
and using ( 16 ), 

((w, (X-L^w)) < sup {2«ti;, /)) - ((/, (A - S)f))} = ((w, (X-Sy'w)) < C((w, (X-So^w)), (46) 
/ 



which is computable because So acts diagonally on 
Subsection 3.1). In the other direction we can use the bound (16) to obtain 

((w, (A - Lr'w)) > Csup {2((w, /)) - ((/, (A - S )f)) - ((Af, (X - So^Af))} . (47) 

This is still not computable because the operators A map Ai n to a (finite) direct sum ®\k\<tM n +k- 
Recall that So denotes the generator of the symmetric simple exclusion (we do not denote the 
dependence on d) . We note that if / = J2a Java then 

So f(«) = E f(v x ' x+ei ) - M = E ^ E - h) ■ (48) 

x,ei A |{a;,x+ei}nA|=l 

Here A x,a;+e ' is the finite set one obtains by replacing x with x+e^ or x+ei with x in A. In particular, 
A = A x ' x + e * unless |A n {x, x + e^}| = 1. So maps the degree n functions, M. n into itself. If we 
think about the coefficients {/a} as a function on finite particle configuration | A| = n then So is just 
the generator of a collection of n particles moving according to continuous time simple symmetric 
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random walk with the exclusion rule. This is the duality property of the symmetric simple exclusion 
process [22] . 

In one dimension So acts on the dimension reduced form of / the following way: 

Sof(yi, Vn-i) = £ V+/(y) + Vr/(y) (49) 
i=l 

where is the effect of moving the i th particle to the right or left (if there is space): 
Vtf(y) = Hy i >0)(f(...,y l -i + l,y l -l,...)-f(y)) 
Vrf(y) = l( w _ 1 >0)(/(...,y i _i-l,y i + l,...)-/(y)) 

where yo = y n = oo. Note that this is a symmetric random walk in WTT 1 where the possible 



moves are ±ei, =L(e2 — ei), . . . , ±(e n _i — e n _2), ±e n _i and each can happen with rate one (if the 
corresponding move is allowed, i.e. we stay in Z? _1 ). The Dirichlet form in the dimension reduced 
form is given by 

1 n 9 

(f,(S )J)= 2 ££(V+/(y)) • (50) 

V i=l 

In any dimension d, So operates diagonally on the direct sum of the M n , and in each is just a 
Laplace operator on the antisymmetric subspace of {'L d ) n . The variational problem involves the 
operators A which couple the different Ai n . Computing the action of A from M. n to Ai n+ k, k>l, 
it takes appropriate differences of a function and employs them as various boundary values for 
the part of the variational problem in M n +k- Thus the problem can be represented as a set of 
variational problems on positive quadrants of (Z d ) n , with the different ns coupled through boundary 
conditions. Even a cutoff version of the problem, optimising over test functions of degree two or 
three, turns out to be challenging. 



If one restricts the supremum in (45) to functions of degree less than or equal to n, one of course 



obtains an increasing sequence of lower bounds. On the other hand, the last term ((Af, (X—S) 1 Af)) 



in ( 45 ) is itself given by a variational formula, and one can restrict the test functions there to obtain 
a decreasing sequence of upper bounds. In the case of the exclusion process at p = 1/2, slightly 
tweaking this reasoning by simply solving an approximate resolvent equation restricted to functions 
of degree less than or equal to n one can obtain an oscillating sequence of upper and lower bounds 



similar to (46) and (47) via a continued fraction- like approximation of (A— L)~ l (see |19] for details). 
This relies on the fact that the generator in this case is tridiagonal on the graded space and that the 
diagonal part is symmetric. For the general speed change generators the upper and lower bounds 
no longer hold. So these estimates are not available in our context. In addition, in the two special 
cases d = 1, j"(0) / and d = 1, j"(0) = 0, the sequence of upper and lower bounds seem to be 
very slowly convergent. The strange case is d = 2, j" (0) ^ which has logarithmically diverging 
diffusivity, but fast convergence of these bounds. In a sense, the case d = 1, j"(0) = is the hardest 
of all, because of the slow convergence coupled with the lack of a solvable model to compare to. 

3. Key tools and proofs of the main results 



3.1. Estimates on the symmetric part of the generator. For the ci-dimensional SSEP we 



have Sq = J2f=i Z)f=i(^ 7 j,-e,- + ^i,ej) where Vi,± ei is the effect of moving the i th particle in the 
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direction ±ej (if possible). The formula for (/, (— So)/) (i n the dimension reduced form) is similar 
to (50): we just have to sum (Vi,± ei f) 2 on B n [d] and i = l,...,d. 

Lemma 3. Let / : B n [d] — )■ M be a fixed function and A E S n [^]- Then 

f(A) 2 <C n (X)(f,(-S )f) (51) 
where there is a C = C(d) < oo such that 

C if n > 4 and d = 1, or if n > 3 and <i = 2 or if n > 2, d > 3, 

C„(A) = I C|logA| ifra = 3,d = l orra = 2,d = 2, (52) 

if n = 2,d = l. 

Proof. We need to consider the dimensions case by case: 

d = 1. As we have already discussed, we may represent B n [l] as TTC . In order to show (51) for 
n > 4, d = 1 we may assume that n = 4. Indeed, we may consider the function / : Z 3 ^ — > M with 
f{x\, X2, x 3 ) := /(xi, x 2 , x 3 , 2/4, . . . , y n _i) for some fixed 7/4, ... , y n -l and apply the n = 4 case for 
/ and x = (xi,x 2 ,x 3 ). The only thing we have to check is that (/, (— S),/) can be bounded by 
C(f, (-S)f) with a fixed C. But this is clear from ( J50| ) and 

V|/(X1,X2,X 3 ) = -Vf /(xi + 1, x 2 , x 3 , . . . ) - /(xi, x 2 + 1, x 3 , . . . ) - /(xi, x 2 , x 3 + 1, . . . ), 

(53) 

which just states that the we can increase the third gap among n particles by moving the first three 
to the left by one step. 

To prove (51) we need to show that 

fl J^2 </.(-W>>0. (54) 

where / : Z+ — > K, y G Let / be the extension of / to Z 3 by reflections: /(dbyi, ±y 2 , ±2/3) := 
2/2) 2/3)- Then 8(/, (—So)/) > (/, (-S)f) where S is the generator of the same random walk, 
but now on Z 3 . The new variational problem can be solved exactly by Fourier transform: 

inf -r\2(9i(-S)g) > C[ { , 2 , 77 7^77: 2 2 dtidt 2 dtA > 0. (55) 

9 £7(2/) 2 V/[-T,T] 3t i + (*2-ii) 2 + (i3-i2) 2 + i 3 / 

For the n = 3 the same argument gives 

M-^{g,{X-S)g)>c( f Ty^ dtldt2 ) >C\log\\-\ (56) 

9 9(2/) 2 V-A-m 2 A + t i + (*2 - hr + t% j 

for the n = 2 case one gets the bound 

iriJ—(g,(\-S)g)>c([ YT^ dt ) >CA^/ 2 . (57) 

9 g[y) \/[-t,t] 2 A + r y 

d = 2. As in the previous case, we may assume that n = 3. We may also assume that A = 
{xi,x 2 ,x 3 } has three elements with different x coordinates (which are ordered according to the 
index). Otherwise we could move two of the particles by one unit, and we could bound the 
difference with a constant times (/, (— So)/). Now consider a new generator S where x\ and 
x 2 performs symmetric random walks with the constraint that the order of the x-coordinates of 
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the three particles is preserved. Clearly (/(—So)/) > (f(—S)f) and S corresponds to a symmetric 
simple random walk on (Z + x Z) 2 . The proof now can be finished similarly as in the n = 1 case. 

For n = 2 the proof is similar, we may consider the generator of the simple symmetric random 
walk on Z + x 7L instead of S and after the reflection trick (see the d = 1 proof) the variational 
problem can be computed explicitly. 

d > 3. The proof goes along the lines of the previous arguments: we may assume n = 2 and 
replace S with the generator of the simple symmetric random walk on Z+ x Z d_1 . □ 

Lemma 4. 1. (H-\ bounds for homogeneous polynomials.) 



d=l 

d = 2 
d > 3 



7a,(A-5o)-V ?a )) < 

7a, (a - s rv ?A » < C 





if 


|A| = 2, 


C\ logA| 


if 


|A| =3, 


c 


if 


|A| >3; 


C\ logA| 


if 


|A| =2, 


C 


if 


|A| > 2 



(58) 

(59) 
(60) 



2. (Triviality criterion) For all d, k > 1, let c\, A C Z rf , |A| = k be a finite collection of coefficients 
with J2a°a = 0. Then there is a C depending on the ca such that 



«£ c A 7? A , (A - So)' 1 E c A f] A }} < C. 

A A 
Proof. By the variational formula for any given a > we have 

((fa, (A - SorV)) = sup 2a((f) A , /)) - a 2 ((/, (A - S )f)). 

f 



(61) 



(62) 



Thus if we have a uniform bound ((r/A, /)) < C((f, (—So)f)) then {{fj\, (A — So) t)a)) < 2Ca — a 
which gives an upper bound = C 2 with the choice a = C . From this most of the statements follow 
from the previous lemma, we only need to prove (61). 

If J^A C A = and all subsets A are the same size then we can represent / = J^A ca^a as a finite 
linear combination of gradients t)q — 7/5 where |0| = |0| = k and SI = Vt x ' y where x, y are nearest 
neighbors. (This is clear if / is of the form t)a 1 — fjA 2 an d the general case follows by linearity.) We 
will show that 

((Wi ~ fin, (A - Sbr 1 ^ - r)n))) < !> 



from this the (61) follows immediately. Using the dimension reduced variables in Bk[d] and the 
variational formula we get 



((vn - f)n, (A - S ) 1 {fjn - %))) = sup2(g A - g~ h ) ~ (9, (A - S )g) 



(63) 



But (g, (A — So)g) > (g\ — 9\) 2 so the right hand side is bounded, which completes the proof. □ 
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3.2. Reduction of the asymmetric part of the generator. For the lower bound computations 
we will only need the evaluate the asymmetric part of the generator on degree two or three functions 
(before dimension reduction). In this section we will show how we can simplify these operators in 
the general case. We will work with the generator in the dual representation, and in the one 
dimensional cases we also use dimension reduction. 

We review the simplest case of the one dimensional TASEP [19] here, before going on to discuss 
the general speed change model. Although the process is totally asymmetric, the generator L is 



the sum of a symmetric part So = k(L + L*) (49) and an asymmetric part Atasep = \{L — L 



The latter can be expressed as a sum A$ + A + + where A + maps a degree n function into a 
degree n + 1 function, Aq = —Aq preserves the degree and A_ = —A* + decreases the degree by one. 
They are written in duality after dimension reduction as 



1 2 n 

(A f)(yi, ...,y n ) =— Hw-i > °)(/(- • • ' W-i - 1, % + 1, • • • ; 
z i=i 



f(y)) 



i 



o n 
2p 



8=1 



i(vi > o)(/(. . . , m-i + - 1, . . . ) - f(v))\ 



(64) 



(A+f)(yi, ...,y n ) = y/x^HVi = 0)(/(. . . ,yi-i + l.j/i+i, ...)-/(.. .,y i ^ 1 ,y i+1 + 1, . . . )); (65) 



8=0 

n 



. . . , y n ) = HVi > 0)(/(- • • , Vi-l + 1, 0, y 4 , ...)-/(... , y^ lt 0, Vi -1,...)). (66) 



8=1 



In all of these formulas yo = y n +i = oo. 

Our lower bounds will rely on finite degree test function computations. More specifically, we will 



use the bound (47) and restrict the supremum to certain finite degree test functions. Recall the 



definition of K from the assumptions on our rate function r( 
In case of the general one dimensional model we will use 



10K 



MP := j/eM 2 : /» = f(y) if 
For the one dimensional model at inflection we consider 

A-4 1} := {/ G M 3 : f{x, y) = f(y, x) for x, y G Z + and 

Vi 



10K 



for x, y £ 



(67) 
(68) 



f(xi,x 2 ) = f{yi,yz) if 



WK 



for i = l,2,xi,yi G Z+,f(x,y) = f(y,x) 



Finally, for the general two dimensional model we will work on 



M { ? := I f G M 2 : £ f A+x = £ f {0 ,l} +x , for every |A| = 2, A C [0, WK] 2 } . (G9) 

Lemma [5] below describes how A simplifies in the general one dimensional cases if we restrict it to 
the specific classes of test functions. It shows that A can be approximated with a constant multiple 



of A + + A^ from the TASEP generator (see ( 65 ) and ( 66 ) ) and the linear combination of certain 



simple operators. The key point is that at the inflection point the TASEP part of the generator 
vanishes. 
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In order to deal with the negligible terms we introduce the following definition. We say that an 
operator T is sectorial on Q if there is a C < oo such that for any / 6 Q, 

((Tf, (A - So^Tf)) < C((f, (A - So)- 1 /)). (70) 

The idea is that any sectorial part of the operator A can be dropped from the right hand side of 



(471) without affecting the quality of the lower bound. Our dream was that A + 3 -4f^- Atasep would 



be sectorial 5 ! This is unfortunately false. What is true is that one can find a reasonably simple 
operator A so that A + - ^ Atasep — A satisfies a graded version of the sector condition, i.e. 
it is sectorial when restricted to a finite degree subspace. Since the bounds one obtains improve 
very slowly with degree, one might as well restrict to the lowest degree where they give something 
non-trivial, which turns out to be degree 3 in d = 1. The operator A is not unique, and we only 
really care about how it acts on / £ .M 2 U M 3 . We will prove that the following choice will 
work: In dimension reduced representation, if / 6 M 2 then 

Af( yi ,y 2 ,y 3 ) =C(l {yi=y2=oy (f(y 3 + 1) - f(y 3 )) + l {w=w=0} (/(ift + 1) - f(Vi + 2))); (71) 

If / G we have A = c 1 A 1 + c 2 A 2 + c 3 A 3 with 

A 1 f(yuy 2 ) = H {yi =o}(/(0,y 2 ) - f(0,y 2 + 1)) + l {y2=0} (f(yi, 0) - f( yi + 1,0)) (72) 
A 2 f( Vl ,y 2 , y 3 ) = l {yi=y2=0} (f(0, y 3 + l)- /(0, y 3 )) + l {w=w=0} (/(yi, 0) - /(yi + 1, 0)) (73) 
A 3 f(yi,y 2 ,y 3 ,y4) = ^{ yi = y2 =o}(f{y3 + 1,2/4) - /(j/3,2/4)) + l{y 2 =y 3 =o}(/(2/i + 1,2/4 + 1) - /(yi + 2,y 4 )) 
+ l{ w=w =o} (/ (2/1, 2/2 + 1) - /(2/i,y 2 + 2)) (74) 

The following lemmas are proved in Section |4j 

Lemma 5 (Simplified form of the asymmetric part in Id). Consider a one dimensional coercive 
lattice gas with local rates satisfying the divergence condition, and denote the asymmetric part by 
A. Let Atasep be the asymmetric part of the generator of TASEP for p = 1/2. Then there is an 



operator A satisfying (71) - (74) so that 



A + I^-Atasep ~ A (75) 



is sectorial on Ai^ U .Mc 1 ^ 



Lemma 6. Suppose that we have a lattice gas in two dimensions satisfying our conditions. Then 

(2) 

A — Atasep is sectorial on Ai 2 where Atasep is the asymmetric part of the generator for a 
nearest neighbor 2d TASEP with drift —2(j'((p),j 2 (p)). 



Note that it is easy to construct a nearest neighbor TASEP with a given drift (a, b). Assuming 
a, b > one can just take the jump law to be p{e\) = a + l,p{—e{) = l,p(e 2 ) = b + l,p(—e 2 ) = 1 
with p(y) = otherwise. 



5 Recall the for TASEP j"{p) = 



—2, hence this is really a difference 
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3.3. Estimates allowing removal of the hard core. Although the previous lemmas show that 
A can be simplified on -M^ , .M3 and even the simplified operators are too complicated to 

do actual computations. The main reason is that because of the exclusion condition (the 'hard core 
interaction') the test functions do not live on Z d , but on Z rf minus diagonals. It is natural to try 
to remove the hard core interaction, i.e. approximate the test functions, and operators, with ones 
on Z d . 

This was one of the main tools in |19| and [37] . If one considers a test function in M n in dimension 
d then one can naturally extend this to a symmetric n-variable function on Z d which vanishes when 
two coordinates coincide. The generators extend naturally (e.g. Sq will correspond to the lattice 
Laplacian) and one can use Fourier techniques to do computations with the new versions of these 
operators. The 'removal of the hard core' lemmas of [19] (see Section 4 in that paper) show that 
when we make these approximations we do not make a large error when we replace ((/, (A — So)/)) 
and ((AtasepI, (A — So)~ 1 ATASEpf}} with their approximate versions. If d = 2 then the extra 
constants appearing in the estimates will only depend on n (and not A) while in the d = 1 case an 
extra | log A| multiplier appears. 

We are still able to use the d = 2 version for our purposes, but the | log A| factor would kill 
our estimates in the d = 1 inflection case, because the order of the diffusivity is smaller, so we 
use a different approach. We consider the dimension reduced form of our test functions to get an 
extension on Z (in case of JVl^ ) and Z 2 (in case of and show that the approximate versions 

of the appropriate quadratic forms give good enough bounds. These are the contents of Lemma [7j 
and Lemma [8] below, whose proofs are given in Section [5] 

Lemma 7. Consider a one dimensional lattice gas model satisfying our conditions. If f : Z + — > R 
is a test function with f € then 

(f,(X-S)f)<C [ (X + t 2 )\F(t)\ 2 dt (76) 

(ATASEpf, (A - S^ATASEpf) < C [ ^=\F(t)\ 2 dt (77) 

J\— 7r,7r] V A 



[— 7r,7r] V A 

{Af, (A - Sr'Af) < C [ (A + t 2 + f-)\F{t)\ 2 dt. (78) 

J[— 7r,7r] vA 

Here A is the operator from Lemma^ and F(t) = J^'k—o f {k)(e ttk + e lt ^~ k ^) is the Fourier transform 
of the following symmetric extension of f: 

F(x) = f(\x\ - l x<0 ) (79) 

Lemma 8. Consider a one dimensional lattice gas model satisfying our conditions. If f : Z 2 ^ — > R 
is a test function with f G then the following bounds hold: 



(f,(X-S)f)<C (X + ti + t^Fih^dhdh (80) 

JJ[-tt,tt] 2 

(Af,(\-S)- l Af)<C [[ (tj + tDllogXWFitut^dhdtz (81) 



+ C 



[— 7T,7r] 



F(t 1 ,t 2 )dt 2 



\h\dti + C{f, (X-S)f). 
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Here A is the operator from Lemma \^ and F(ti,t 2 ) is the Fourier transform of the following 
symmetric extension of f: 

F(x, y) = f{\x\ - l x<0 , \y\ - 1„< ) (82) 

Note the previous lemmas will be used for dimension reduced functions, and therefore are stated 
with only the single bracket. 

Lemma 9. Consider a two dimensional finite range exclusion model with drift ae\ + be 2 - Let f be 
a test function in M.^ and assume that F : 1? — > M. is an extension of f, 

E /({fa + W,S2 + ya), (2/1,2/2)}), if (xi,x 2 ) 7^(0,0), (83) 
2/1,3/262 



F{Xl,X2) 

with F(x, x) 



F(-Xl, -X2) 

0. Then 



[— 7r,7r]- 



(\ + 4 + tl)\\F(t 1 ,t 2 )\ 2 dt 1 dte 



((f, (A-S„)/»<C 

+ A_)/, (A - S rV+ + 4_)/» <C II (ah + 6t 2 ) 2 | log A||F(ti, t 2 )| 2 dti<ft 2 . 

JJ[-7T,7r] 2 

i^ere ^4+ + j4_ is i/te asymmetric part of the generator of the exclusion without the Aq part. 
3.4. Proofs of the main results. 



(84) 
(85) 



Proof of the upper bounds in Theorem [7} We begin by proving ( 30 ) . By ( 46 ) it is enough to prove 
((vi, (A — So)~ 1 Vi)} < C. By the definition of V{ and W{ it is clear that if we write Uj = J2a c a^a- 
then for any fixed k we have J2\A\=k C A = 0. Hence (30) follows from (61). The rest of the upper 
bounds follow immediately from ( |43[ ) , ( |44| ) and the two parts of Lemma|4 □ 

Proof of the lower bound in Theorem^ in d = 1 when j"(p) ^ 0. We will evaluate the right hand 
side of (47) with an appropriately chosen / G M 2 . By the Cauchy-Schwarz inequality and Lemma 
[5] for any / E JA^p we can bound {{AS, (A — So)~ 1 Af}} above by a constant multiple of 

{{ArASEpf, (A - S Q )- x A T A8Epf)) + Uf, (A - So^Af)) + ((/, (A - S )f)) (86) 



where A is described in (71)-(74). (Since we will work on M 2 we w * n only need (71).) Using 
the dimension reduced variables we may represent / as a Z + — > M function / which is constant on 
intervals of the form [iK, (i + \)K — 1] where K = 1QK. Then by (86) and Lemma M we have 



((/, (A - So)/)) + {{Af, (A - So)' 1 AS)) <C (A + \-W#)\F{t)\ 2 dt 



(87) 



where F(x) is the symmetric extension of /. (We can assume that A is small enough.) 
Since / G Ai 2 we have 

{{w,S))= E "^/(x), where w = ^wa']a and w x = E W A- 

xeZ+ A A={y,x+i/} 

Because of the local evolution condition on r(-, •) and the definition of w we have w x 
x > K = 10-fC and since / G A<2 * ms means that 



««>,/» = 7(0) E ^ 



/(o) E 

|A|=2 



X 



j"(p)7(o) 



for 



(89) 
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where the last identity follows from (43). 
Now choose F : Z — 



so that it's Fourier transform is given b; 
k-i 



G(t) 



1 



l(\t\ < 1/2) 



(90) 



F(t) = d(2Kt) J2 eijt > 

j=-K 

It is not hard to see that F(x) = F(—x — 1) and that F is piecewise constant on intervals of the 
form [jK, (j + 1)K — 1]. Moreover, 



no) 



F(t)dt > C R X-^ 4 



(X + X-^H^Fit^dtKC^X- 1 ^. 



(91) 
(92) 



We assumed j"{p) > in the first line, but we can do that without the loss of generality. (Otherwise 
use — / for the test function.) Now if we choose our test function / to be the restriction of aF(x) 
to x > then by the previous estimates (and p7|)) we have 

((w, (A - L)-^)) > CiaA- 1 / 4 - C 2 a 2 X- 1/4 (93) 

and the desired lower bound follows by choosing a = C\j2C 2 . O 

Proof of the lower bound in Theorem^ in d = 1 when j"(p) = 0,j"'(p) ^ 0. By (42) we have 
J2\A\=2 W A = but X^|A|=3 7^ 0. By (58) and the Cauchy-Schwarz inequality, if we let w = 
J2\M=3 C ^VA, we have ((it?, (A — L) _1 u;)) > \{{w, (A — L)~ l w)) — C. We use the variational inequality 
(47) with w instead of w and get the lower bound by choosing an appropriate test function / £ .M3 . 
Using the dimension reduction (see (38)) we may represent / as a Z^_ — > R function. By the Cauchy- 
Schwarz inequality and Lemma [HJ we can bound ((Af,(X — So)~ 1 Af)) by a constant multiple of 
((Af, (A — So)~ 1 Af)) + ((/, (A — So)f)), since the Atasep part vanishes. By Lemma [8] we have the 
upper bounds 



((/, (A - S )f)) + ((Af, (A - Sor'Af)) < C 

+ C 

+ C 



2 (X + t 2 1 + t 2 2 )\F(t 1 ,t 2 )\ 2 dt 1 dt 2 
(t? + i2)|logA||F(ti,t 2 )| 2 dtidt 2 



(94) 



[ — 7T,7r] 



F(t 1 ,t 2 )dt 2 



\h\dh. 



where F is the symmetric extension of / to Z \ described in ( 82 ) . Now choose F so that its Fourier 
transform is given by 



K 



F(h, t 2 ) = G(2Kt u 2Kt 2 ) J2 expiiktx + ijt 2 ), 

k,j=-K-l 

G(h,t 2 ) = 1(A < tj + t 2 2 < 1/2)- _J 



{t( + 1%) \og(t( + 4) 



(95) 
(96) 



°This is similar to the test functions used in [19] 
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with K = 10-KT. Clearly F has the desired symmetry given by (82), it is constant on the regions 
[xK, (x + 1)K — 1] x [yK, (y + 1)K — 1], x, y £ Z + and also F(x, y) = F(y, x). Thus if we choose / 
to be the restriction of F to then / £ .M3 and we have the upper bound (94). One can now 
explicitly bound the integrals on the right side of ( |94[ ) to get 

((/, (A - S )f)) + ((Af, (A - Sor'Af)) < d log | log A|. (97) 

Using the same argument as the one leading to (|89|) we have 

F(tx,t 2 )dtidt2 > C 2 log | log A|. (98) 



«t5,/» = /(o,o) £ ^A = ^/'(p)/«).o)-r 



|A|=3 



3! 



(We assumed j"'(p) > 0, otherwise use — / as the test function.) Choosing af as the test function 
in ( [47] ) with a = C 2 /2C\ completes the proof. □ 

Proof of the lower bound in TheoremYn in d = 2 when j"(p) / 0. We use the same strategy as be- 
fore, here we will use a test function from Ai 2 in the inequality (47) to get a lower bound. Suppose 
that F : I? -)• E satisfies 

F(a;i,2!2) =-F(-xi,-x 2 ), and F(xi, x 2 ) = F(0, 0) if |x 2 | < = 10K. (99) 

Define the test function / € _M 2 as 

F(xi,x 2 ) if (3/1,2/2) = 0, 

otherwise 



/({(xi,x 2 ), (yi,y 2 )}) 

it is easy to check that / G A^^. Then by Lemma ^ and Lemma [9] we have 

(X + tj + tl)\\F(ti,t 2 )\ 2 dtidt 2 



(100) 



((/, (A - So)/)) + «4f, (A - ^o)- 1 ^/)) < C 

+ C 

< a 



[— 7T,7r] 2 

(oti + bt 2 f\ log A||F(ti,t 2 )| 2 dtidi 2 



(101) 
(102) 



[ — 7T,7r]' 



(A + (bh - at 2 y + I log X\(ah + b 2 t 2 ) 2 )\ \F(h, t 2 )\ 2 dtidt 2 



Here (a,b) = (j"(p), j 2 (p)) and C does not depend on A. If we choose 



K 



F(t 1 ,t 2 ) = G({2K + l)t 1 ,(2K + l)t 2 ) J2 expiiktt + ijh), 

k,j=-K 

G{h,t 2 ) = l(tl + t 2 2 < 1/2) (A + (bh - at 2 f + I log A|(aii + b 2 f 



(103) 



(104) 



(105) 



then F(xi, x 2 ) = F(— xi, —x 2 ) and F( x±, x 2 ) = F(0, 0) if |xi|, \x 2 \ < K. We can explicitly estimate 



the integral on the right hand side of (103) and show that it is bounded by C\^J\ log A|. Using the 
same argument as the one leading to (89) we also get 



((m,f)) = ^jUp)F(o,o) = ^f{( P ) 



[[ Ffat^dhdh >C 2V ^oiA. 



(106) 



The statement now follows by choosing af as the test function with a = C 2 /2C±. Note that we 
assumed that j'{(p) > 0, but we may do that without the loss of generality. □ 
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4. Proof of the reduction lemma 



Using Q and A = (L + L*)/2 we have 

^/ = JVxEE^^^+i/ - Vx)(f(r) x ' x+V ) ~ m) (107) 

x y 

Let / = J^aIaVA be local. We can also express r(y, 77) in our basis as r(y,rj) = X)a c ?/,aOa, and 
because the rates are assumed to be local, each A is a subset of [-K, K) d and An {0, y} = 0. Using 
the identity 

fg = l + Kf) x1 K = (l-2p) X ' 1/2 , (108) 

we then have A = J2y,A ^A,y where 

M,yf = lVXCy,A Y Y VA+x(2 + nfj x+y + K7) x - 2<f){ XjX+y y)(f niJ { x } - fnu{x+y})vn- (109) 
x nn{x,x+y}=(A 



By (34) and (108) if and are disjoint then fjofjn = r/eun, otherwise we have rye 5 ?^ = VBAn + 
Kfjeun where XAY denotes the symmetric difference. Hence (109) becomes 



2 VXCy,A J2) Y ( 2 + Kf lx+y + KVx ~ ^V{x,x+y} ) (fi 

x I fin(A+x1=0 



fiu{i} — fnu{x+ y })vnu(A+x) (HO) 



nn(A+x)=Hi 

Qn{x,x+y}=% 



+ Y Y ( 2 + K ^x+y + Kfj x - 2f){ XjX+y y)(fn u { x } - /sj U {i+y})( J ?siu(A-n)\(r+i) + K Vnu(A+x)) 

0^rcA nn(A+x)=r+a; 
nn{x,x+y}=0 

(111) 

It becomes clear one needs some simplifying notation. If B\ and B% are disjoint finite subsets of 
Z d with 0, y £ B 1 U B 2 and B 3 is a subset of B\ U i?2 with y £ B%, ^ -B3 we let 

A[5i,S 2 ,B 3 ,y]/ n = Y ^{B 1 +xCQ,B 2 +xn(l=9} (/(n\(B 1+a: ))ij(B°' I '+s) ~ /(n\(Bi+x))U(B 3 +x)) (H 2 ) 

where B + x means that we add x to every element of B. In plain words: 

1. For each x £ Z d we check if the shift of Q matches the pattern given by B± and B2 in the sense 
that Q contains B\ + x and is disjoint with B2 + x. 

2. If it matches for a given x, then we replace the part B\ + x in f2 by B% + x and B 3 ,2/ + cc 
respectively, evaluate / at these two values and compute the difference (B 3 ' v means that we delete 
y and add the element to B3) 

3. We sum over all the x where we have a match. 

If 17 = {p\, . . . ,p n } in lexicographic order and the smallest element of B\ is q. Then we can 
rewrite the previous line as 



Y 1 {B 1 +(p i -q)cfl,B 2 +(p i -q)rfo=<i} (/(n\( J Bi+ a ;))U(B°' ! '+(p. 1 - g )) ~ /(fi\(Bi+(Pi-q)))U(B 3 +(p i -g)) ) (H 3 ) 



1=1 
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We can also compute the dimension reduced form of this. Let A £ B n [d] and assume that f2 
{pi, . . . ,p n } is in the equivalence class A. A simple computation gives 



A{B 1 ,B 2 ,B 3 ,y]f A - l{Bi+(p»-<z)cfi and B 2 +( Pl -g)nn=0}(/A l - /aJ- 



(114) 



i=l 



Here / is the dimension reduced form of / (see (38)), and Aj, Aj are the equivalence classes of the 
sets (Q \ (Bx + x)) U (B®' y + (p» - q)) and (Q \ (B 1 + (p* - q))) U (S 3 + (p; - g)), respectively. Note 
that this can now be written as the sum of operators of the form -A[i?x, B 2 , B3, y]. Indeed, the sum 
(110) is just \-\fx°y,A times 

2A[A, {0, y}, {y},y} + kA[A U {y}, {0}, {y}, y] + kA[A U {0}, {y}, {y} 5 y] - 2A[A U {0, y}, 0, {y}, y]. 

(115) 



The sum ( 111 ) will have similar form, for each / T C A we get \^fx c y,A times the following sum: 
2kA[A, {0, y}, T U {y}, y] + k 2 A[A U {y}, {0}, T U {y}, y] (116) 
+ K 2 .4[A U {0}, {y}, T U {y}, y] - 2kA[A U {0, y}, 0, T U {y}, y] 

2A[A \ r, r u {0, y}, r u { y }, y ] + «a[a u {y} \ r, r u {o}, r u { y }, y ] (117) 
+ kA[a u {0} \ r, r u {y}, r u { y }, y ] - 2A[A u {0, y} \ r, r, r u {y}, y ]. 

We will show that many of the operators A[B\, B 2 , B^,y] can be bounded using the Dirichlet form, 
and that in a certain sense the important quantities about such an operator are \B\\ and I-B3I. The 
following lemma summarizes what we have proved so far. 



Lemma 10. Consider a fixed finite A with \ A\ = I andy 6 Z . Then the operator A ^ y of (109) can 
be expressed as a finite linear combination of operators of the form A^ y = otB!,B2,B 3 A[Bi, B 2 , B^,y]. 
The sum of the coefficients aB 1 ,B 2 ,B 3 with \Bi\, \B 2 \, \B$\ fixed is given by the following table: 



\Bi\ 


\B 2 \ 


\B 3 \ 


sum of coefficients 


£ 


2 


1 




£ + 1 


1 


1 




£ + 2 





1 


~VXCy,A 


£ 


2 


o + l 




£ + 1 


1 


o + l 


^y/XCy^O 


£ + 2 





a + 1 




£-a 


a + 2 


a + 1 


K VXCy, aQ 


l + l-a 


a + 1 


a + 1 




£ + 2-a 


a 


a + 1 





;ii8) 



where 1 < a < I. 

Examples. 1. Exclusion process. We consider the Id TASEP first, in that case the only non-zero 
Cy : \ is C10 = 1. Then A + ,A^,Aq can be expressed as 

K 



A 
A + 



X(A[{1},{0},{1},1]+^[{0},{1}, {!},!]) 



XA[{O,1},0, {1},1] 



A- 



XA[0, {0,1}, {!},!]. 



(119) 
(120) 
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In case of a general exclusion with jump law p(-) the operator A is just a linear combination of 
terms like this with y in place of 1 if p(y) > 0. 

2. The model from (18). For the model (18) at density p = 1/2, k = 0, so many terms vanish. 
The rates in our basis are r(l,rj) = 2 — hfi-i — \r\i-, r(— l,rj) = 2, and only the — |f)_i — J772 part 
contributes to the asymmetric part. This means that 

4A =A[{-1}, {0, 1}, {1}, 1] - A[{-1, 0, 1}, 0, {1}, 1] + A{{2}, {0, 1}, {1}, 1] - A[{0, 1, 2}, 0, {1}, 1] 
+ A[0, {-1, 0, 1}, {-1, 1}, 1] - A[{0, 1}, {-1}, {-1, 1}, 1] (121) 
+ A[0, {0, 1, 2}, {1, 2}, 1] - A[{0, 1}, {2}, {1, 2}, 1] 



Lemma 11. Suppose B\,B 2 , B 3 , y and B[, B' 2 ,B 3 ,y are as in the preamble to (112) and BiL)B 2 = 
B[UB 2 = [-£,£] d . If\Bi\ = \B[\ then A[B 1 ,B 2 ,B 3 ,y] - A[B[, B' 2 , B 3 ,y] is sectorial. 

Proof. We will prove sectoriality first with respect to the scalar product (•, •). 

We have to bound (/, (A[Bi, B 2 , B 3 , y] - A[B[, B' 2 , B 3 ,y])g) by Dirichlet forms. Using the defi- 
nition of the operators and rearranging the corresponding sums, it can be written as 



£ £ (/ru(Bi+aO - fru(B[+x))(9 
x rr\([-e,e] d +x)=Q) 



ru(B°' y +x) 



9ru(B 3 +x)) 



122) 



Since \Bi\ = \B[\ we can get from B\ to B[ by finitely many nearest neighbor steps (all of which 
take place inside [—£, £] d ): £1® = B\, Q\, . . . , fi™ 1 = B[ where we get f^ +1 from £1\ by moving one 



of its elements by a unit vector. Thus we may rewrite fru(B 1 +x) 
gradients: 

mi — 1 



fru(B[ 



+x) 



fruBx - /rus; — z2 f] 



ru(n\+x) fru(Q\ +1 +x) 



i=Q 



as a sum of discrete 



(123) 



Since \B 3 



\B 3 ' V \ we have an analogous representation for 5 f r u(B 0,!/ +x) — 9ru(B 3 +x) with a series 

B 3 ,nj--- ^r = Bj' y 3 ~ 

thus bounded above by ^fh{f)D^g) where 



of sets ^2 



B 3 ,y By Cauchy-Schwartz inequality the left hand side of (|122|) is 



mi — l 

x rn[-£,e] d +x=9 i=o 



ru(n[+x) fru(n\ +1 +x) 



-x)> 



(124) 



and D 2 {g) is defined similarly. The right hand side of (124) can be bounded by constant multiple 
of (/, (—So)f) (with a constant depending only on £) since each term is of the form (/a — f'j^) 2 
(where A' is obtained via a nearest neighbor step from A) can only appear finitely many times in 
that sum. This is because the bond where the jump happened between A and A' must be inside 
the cube [-£,£] d + x. The same bound holds for D 2 (g) with (g, (-S)g) which yields 

(/, (A[B%,B 2 , B 3 , y] - A[B[,B' 2 , B 3 , y])g) 2 < C e (f, (S)f)(g, (-S)g). (125) 

This shows the sectoriality for the scalar product (•,•). The statement of the lemma will follow 
from (28) which shows how to get the appropriate bound for ((•,•)) from (•,•). □ 



Lemma 12. Assume that B\UB 2 = [—£,£] , and B\, B 2 , B 3 ,y are as in the preamble to (112) 
Then if d = 1 and \Bi\ > A, or if d = 2 and \B 3 \ > 3, then for T = A[Bi, B 2 , B 3 , z] and f G M r 
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we have 

((Tf, (A - Sy'Tf)) < Cn((f, (-£)/)) (126) 

where C only depends on £. 



Proof. We will use the notation hi = \B{\. From the definition ( 112 ) it follows that if / 6 M n , n > 63 
then Tf E Ma with n = n — 63 + h\. We can clearly assume that n > 63, n > 61. 

We will start with the proof of the one dimensional case . By Lemma [TT] we may assume that 
B\ = [—£, — £ + 61 — 2] U {^}. In the dimension reduced form our operator T = A[B\, i?2, -B3, y] is 

h-bi + l 

I7(3/i, • • • ,j/n-i) = J] = = • • • = ?/i+fei-3 = 0,y i+;)1 _2 = 2i-k)(f(jly) - f{ny)). 

i=i 

(127) 

Here r^y and r ? 'y are defined the following way. Let cZbe the size h set whose dimension reduced 
form is y = (3/1, . . . , j/n-l) an d its i*' 1 largest element is — £ (this will uniquely determine O). Then 
Tjt/ and Tj'y are defined as the dimension reduced forms of (Q \ B{) U B3 and (f2 \ 5i) U -B^' 2 - We 



will use the shortened notation 1, for the indicator function in ( 127| ). By the variational formula 



ft— 61+1 

((Tfi-Soy'Tf)) = sup X 2((3,li(/o^-/oT i )»-((0,(-W> 

Consider the symmetric simple exclusion generator Sf )1 on Mb 1 in the dimension reduced form, 
i.e. acting on functions of the form g : — > E (see (49)). Now let be the operator Sb 1 



acting on dimension reduced variables yi, yi+i, ■ ■ . , Vi+bx-2 of a A^n function. This is the generator 
of symmetric simple exclusion performed by the particles i, i + 1, . . . , i + h\ — 1 where particle i 
is 'glued' to the particles 1, 2, . . . , i — 1 and i + 61 — 1 is glued to the particles i + fei, . . . , h: i.e. 
whenever i or i + 61 — 1 jumps, the corresponding particles perform the same jump simultaneously. 
We claim that for any g S Ai n we have 



«<?, - V S lM g)) < 2n((g, (S )g)) (128) 




The operator Ya=i 1+1 61 corresponds to the dimension reduced form of a random walk where we 
can move each particle left or right if possible with rate at most h\ or we can move the leftmost k 
or rightmost k (1 < k < h — h\ + 1) particles together one unit left or right with rate one. Note 
that whenever such a simultaneous jump happens we can recreate this jump by taking at most 



h — h\ + 1 nearest neighbor jumps with the particles. From this the upper bound (128) follows by 
the Cauchy-Schwarz inequality. 
By Lemma [3] we have 

{{g,-S jM g)) > C((l i5 ,l i5 }}. (129) 

Here we first fix the variables yi, . . . , Vj^i-.yj+bi-iiUj+bx Uh-i, apply the lemma and then av- 
erage out in the variables we fixed. 



DIFFUSIVITY OF LATTICE GASES 



25 



Using this with the bound (128) and the usual variational formula for ((g, (—So) 1 g)) we get 

«T/, (Soy'Tf)) < sup Xf (2((g, lj(f o rj - / o Tj ))) - ~ « 5 , -S^))) (130) 



9 i=l 

n-bi+l 



<sup ^ 2p i5 ,l J (/orj-/or i )))-^((l^,l i5 )) (131) 
9 i=i 

<2C7n £ p i (/orj-/oT i ),l i (/orj-/or i )}) (132) 

<C"n((/,(-5 )/)). (133) 

Here the last step is just Cauchy-Schwarz again: for each j the difference / o rj — / o Tj describes a 
size z jump of a particle which can be recreated by z nearest neighbor jumps. Each such jump can 
appear for at most b\ of the indices j from which the last inequality follows with a C\ depending on 
bi, z. Since n < cn with a suitable c depending only on b%, 63 the statement of the lemma follows. 



Now we will turn to the two dimensional case . Using Lemma 11 we may assume that (—£, —£) 
and (£,£) are elements of B\. We will work with the dimension reduced picture, recall that B n [2] 
is the set of equivalence classes of the size n subsets of Z 2 under shifts (A ~ A + x, x G Z 2 ). 

Let pi denote the i th element of A G Bn[2] in the lexicographic order (this is well defined as the 
order does not change when we shift the set). Then we have 

«T/(-S )- 1 T/» = su P 2((T/, 5 )) - {(g, (-S )g)) 
9 

h—bx+l 

= sup £ 2((g,Ufor' l -foT l )))-((g,(-S )g)) (134) 

3 i=l 

Here li(A) is the indicator of the event that if we consider the shifted version of the box [— £, £] 2 
whose lower right corner is exactly pi then the intersection with A is exactly the appropriate shifted 
version of B\. (This is basically the same as the lj in the one-dimensional case.) The operators 
Ti,r- are just the analogues of their one dimensional counterparts: we replace the intersection of 
the box and A with the shifted versions of B% and 

The rest of the proof is similar to the one dimensional case: we will eventually prove that for 
any g G Mn we have 

n-61+1 

Cn((g,(-S )g))> E Ptf,ltf»- (135) 
i=i 



From this the statement will follow exactly the same way using (134) and the end of the one 
dimensional argument. 

Fix 1 < i < h — 61 + I and suppose that for A G Bn[2] we have Hi (A) = 1. We can move the i — 1 
particles which are in front of pi in the lexicographic order so that their z-coordinates are smaller 
than that of pi. Clearly we can do this by moving each of these particles one step to the left. We 
can also move the particles which are not in the square Pi + [0, 2£] 2 and not in front of pi so that 
their x-coordinates are bigger than 21 plus the x-coordinate of pi . This can be achieved by moving 
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each such particle 2£ + 1 steps to the right. Denote the new configuration by <Tj(A). (If Hi (A) = 
then we may define <7j(A) as the empty set.) 



5, 



Red 



Blue 



Green 



Figure 1 . Schematic representation of &i 

Note that because we made at most 21+1 horizontal steps with each particle to go from A to 
<jj(A) we have 

((Hg -goat), Hg - g o < C(£)((g, (S )g)). (136) 

We will now color the particles of cx^A with red, blue and green so that the particles that are in 
Pi + [0,2-£] 2 initially are red, the particles which have smaller x-coordinates than p, are blue and 
the rest are green. Now define Sibi similarly as before: it will be a symmetric simple exclusion on 
the b\ red particles, but whenever we would change the smallest red x-coordinate then we perform 
the same change on all the blue points simultaneously (i.e. move them one step to the left or right) 
and whenever we would change the largest red x-coordinate then we will perform a similar step on 
all the green particles. It is clear from Lemma [3] that 

((l,9°^li»o4<%(-%)s» (137) 
We also have (by the same argument as in ( 128[ )) 




L S ijbl \g))<2h((g,(-S )g)). (138) 
Putting together our bounds will yield (135) (with an £ dependent constant C) and this completes 



the proof of the lemma. □ 
Remark 13. If ^ z <£ B x U B 2 then 

A[B 1 ,B 2 , B 3 ,y] = A[B U B 2 U {z}, B 3 , y] + A[B t U {z}, B 2 , B 3 U {z},y}. (139) 



This follows immediately from the definition (112), one just has to rewrite the indicator on the 



right hand side as l{(B lU {z})+a;cn and B 2 +xnn=9} + ^{B^xcn and (B 2 u{z})+xnn=(&}- 

Note that this means that we do not have to assume that B\ U B 2 = [— £, £] d Lemma 
suppose that d = 1, B\ U B 2 C [—£,£] and \B\\ > 3. Then by repeatedly applying (139) we 



12 



Indeed, 



may rewrite A[B\, B 2 , B 3 ,y] as a finite linear combination of A[B[, B' 2 , B' 3 ,y] where B\ C B[ and 
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B} U = [ — One can now apply Lemma 12 for each term to get the sectorial bound for 
A[Bx, B2, -B3, y]. A similar proof applies in the d = 2 case. 

Lemma 14. Let T = A[{z}, {0}, {z}, z] with some 0^ z€Z d . Then if f G M n then 

((Tf, (A - Sor'Tf)) < Cn 2 ((f, (A - S„)/». (140) 



Remark 15. Using the representation (119) with Lemmas 11 and 14 we immediately get that the 
TASEP A is sectorial on M n . 



Proof of Lemma 14.. We first deal with the case d = l,z = 1, the general case will be similar. Note 
that on Ai n the operator T is exactly the generator of an exclusion process (with n particles) with 
rate one unit jumps to the right. We can also consider this as the generator of a random walk on 
the elements of i3 n [l]. For an element A 6 B n [l] denote by r^A the effect of moving the k th largest 
particle one step to the right if possible. In the dimension reduced form this can be written as 

/ s J (3/1 > • • • ,Vi-i + - • • • ) if y £ > 

Tk{yi, ■ ■ -,2/n-l) = < . , , (141) 

[{yi, ■ ■ ■ ,y n -l) otherwise 

We denote by the effect of moving the first k particles one step to the left (with o"o being the 
identity). In the dimension reduced form this is just changing y^ to y^ + 1: 

0"fc(yi, • • -,2/n-l) = (2/1,2/2, • • -iVk-uVk + hVk+l, ■■■)■ ( 142 ) 

Note that T(_ag = and in the dimension reduced form a n is the identity. (Note that r and a 
are different from the operators used in the proof of the previous lemma.) By definition we have 

n 

Tf = J2(f°r l -f) (143) 
i=i 

The last term in the sum can be rewritten as 

n— 1 n— 1 

f°T n -f = fo cr„_i - / = ^(/ o cr„_i - / o cr n _i_i) = ^(/o fj n _i - / o r n _j o a n _i) (144) 

i=i i=i 

which gives 

n— 1 n— 1 

T / = E(/ or *-/)-E(/ oT * o ^-/°^) ( 145 ) 

i=l i=l 

Computing ((g, Tf)} we get 

n— 1 n— 1 

n— 1 n— 1 

y «=1 y i=i 

= E E(/N/) - /(y))G/(i/) - g^mivi > o) 

2/ i=l 



DIFFUSIVITY OF LATTICE GASES 



28 



since (f(ny) - f{y))l(yi 



0) = by definition. Note that if yi > then 

(yi, ■ ■ -,yi-i,yi - • • • ) = t it 2 ■■■ny 



and 



(g(y) - gia^y^Hyi > 0) 
which means that 

((9,Tf)) = 



■Mm 



n-1 



0) J2(9{TjTj+i • ■ ■ ny) - g(T j+1 T j+2 ■ ■ ■ ny)) 

3=1 



y i=l j=l 



(146) 



(147) 



(148) 



where y. . 



Tj+iTj+2 ■ ■ ■ ny. Applying the Cauchy-Schwarz inequality and noting that each / 
gradient is multiplied by the sum of i < n g gradients we get 

n—l n— 1 i 

((g,Tf)) 2 < £EUW - f(y)) 2 x £ E^imj) - g(y hJ )f 

y i=l y i=l j=l 

n—l n—l 

< ™ 2 E E(/M - /(=)) 2 x EE(9(i)-s(i)) : 

2/ i=l y i=i 

= 4n 2 ((f,(-S )f))((g,(-S )g)) 



(149) 

(150) 
(151) 



from which the lemma follows for d = 1, z = 1 by taking g = (A — So)~ 1 Tf. 

If d = 1 and z > 1 then basically the same proof works. In this case the exclusion with jumps 
z on Z can be decoupled into exclusions on the sub-lattices zX, z7L + 1, . . . zTL + z — 1 and on each 
of those subdattices we can apply the z = 1 result. Similar trick works for d > 1: we can decouple 
the exclusion on Z d with fixed jumps of size z into exclusions on one-dimensional subdattices of 
the form a + zTL and we can again use the d = 1, z = 1 result. □ 



Lemma 16. 1. Assume that B\, B 2 , B 3 , y are as in the preamble to (112) with d = 1 and we have 
B 1 ,B 2 ,B 3 C [-K,K],\y\ < K. If \B 3 \ > 1 then A[B 1 ,B 2 ,B 3 ,y] will vanish on M 2 l) after the 
dimension reduction, i.e. for any f G M { 2 ] and g we have ((g,A[Bi,B 2 ,B 3 ,y]f)) = 0. Moreover, if 
\B 3 \ > 2 then A[B\, B 2 , B 3 ,y] will vanish on A4 3 after the dimension reduction. 



2. Assume that B\, B 2 , B 3 , y are as in the preamble to (112) with d = 2 and we have B%, B 2 , B 3 C 
[— K, K] 2 , \y\ < K. If \B 3 \ > 1 then A[B±, B 2 , B 3 , y] will vanish on A4^ after the dimension 
reduction. 



Proof. From definition ( 112 ) it is clear that A [B±, B 2 ,B 3 , y] vanishes on M n if n < \B 3 \. In order to 



prove the statements of the lemma we just have to understand Af-Bi, B 2 ,B 3 ,y] on M\b 3 \ after the di- 
mension reduction. If / G -M|B 3 | then Aj and Aj in (114) are always equal to the equivalence classes 



of B®' y and B 3 , respectively. Thus if / always gives the same value on these equivalence classes then 
A[Bi, B 2 , B 3 ,y]f is always zero which will make the scalar product {{g, A[B\, B 2 , B 3 , y]f)) vanish 
as well by (38). From the definitions (67)-(69) oi M 2 ,A4 3 and M 2 the l emma now follows. □ 



Proof of Lemma [5| Write the asymmetric part of the generator A as a finite linear combination 
A = EcB 1 ,B2,B s>y A[B 1 ,B 2 ,B 3 ,y] with B, t C [-K, K] and \y\ < K. 
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16 



By Lemma 12 and Remark 13 all the terms with > 3 are sectorially bounded. By Lemma 
we do not have to worry about the terms with | JB3 1 > 1 on Ai^ and the terms with |2?3| > 2 
on because they will vanish after the dimension reduction. We will now treat the various 

remaining terms according to the size | B3 \ — \B\\. 

Case 1: |Bi| — |Bs| = 2. Since \B\\ < 3 the only possibility here is \B\\ = 3, I-B3I = 1. We will 
show that these terms can be written as a constant multiple of 1, 0, 1}, 0, {1}, 1] and a term 
which is sectorial on M2 U M3. As a first step we note that we may assume that B\ U B2 = [—£, £] 
for some I, one just needs to apply the argument in Remark 13 and use the fact that the building 
blocks A[B[, B' 2 , B 3 ,y] with \B[\ > 3 are all bounded. Next we can assume that y = 1: in the 
definition ( [mj ) we can just replace the gradient 

AfA^i+^MB^+z) ~ /(fAC-Bi+aOM-Bs+af) ( 152 ) 

by a sum of gradients involving nearest neighbor switches B3 — > B^' z 1 . Then by Lemma [ll] we 
may assume that B\ = {—1, 0, 1} so we only need to deal with ^4[{— 1, 0, 1}, B2, {1}, 1]. Finally, by 
using ( 139) repeatedly we can write A[{ — 1, 0, 1}, B2, {1}, 1] as a sum of A[{ — 1, 0, 1}, 0, {1}, 1] and 
finitely many building blocks of the form 1, 0, 1} U T, B2 \ T, {1}, 1] with |T| > which are all 



sectorial by Lemma 12 



From the definitions it is easy to check that ^4[{— 1, 0, 1}, 0, {1}, 1] acting on .A/f 2 and can be 
written as ( 71 ) and ( 74 ) in the dimension reduced form. 

Case 2: |Bi| — | B3 1 = 1. We will need to deal with two sub-cases: |-Bi| = 2,|I?3| = 1 and 
|£>i| = 3, |i?3| = 2. We first start with the case = 2, \B^\ = 1, in that case B3 = {y}. 
Suppose first that B2 7^ 0. Repeated use of ( 139[ ) gives 

A[B 1 ,®,{y},y]= ]T A[B X L)B,B 2 \ B, {y} U B, y] (153) 

BcB 2 

which means that A[B\, B2, {y}, y] — A[B\, 0, {y}, y] can be written as the sum of building blocks 
with = 3, 1-B3I = 2 or \B\\ > 3. Thus that we may assume that B2 = and our operator is 
A[{0, y}, 0, {y}, y\. Note that x 1 ^ 2 ^4[{0, y}, 0, {y}, y] is exactly the A + part for an exclusion with 
jumps of size y (with rate 1). By Theorem 5.2 of [35], the generator of any local mean zero exclusion 
process is sectorial. Applying this with / £ -M n ,g G M-n+i we get that the previous inequality 
also holds for just the A + part of the generator as well (and then it also holds for the scalar 
product ((-, •))). Since vl[{0, y}, 0, {y}, y] — yA +j TASEP is the A + part of the generator of a mean 
zero exclusion (up to a constant multiplier) this implies that we may replace A[{0, y}, 0, {y}, y] with 
yA +t TASEP- This shows that the contribution of the terms with |£>i| = 2, |i?3| = 1 can be written 
as a sum of c+A + with some constant c+, A[B±,B2, B3, y] with |£>i| = 3, 1-B3I = 2 and a term which 
is sectorial. 

In order to determine the value of c+ we just need to collect all the coefficients of the = 
2, {B^l = 1 terms (weighted by x _1//2 2/)- We do this by consulting the table (118). Note that 
|i?i| = 2, |i?3| = 1 can happen if I = 0, 1 or 2 and we get 

c + = y c v,® ~ K v c vA- y°y, A - ( 154 ) 

y y,\M=i y,\A\=2 

To summarize: the contribution of all the terms from Case 2 can be written as c+^4+, ^4[-Bi, B2, B3, y] 
with I Si I = 3, IB3I = 2 and sectorial terms. 
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Now we turn to the \B\\ = 3 
dimension reduction by Lemma 



16 



2 terms. Note that on Ai^ these will vanish after the 



B 3 \ 

To simplify the \Bx\ = 3, 1-B3I = 2 terms onA^j ' we may assume 
that y = 1 (by the same argument as in Case 1). Then B 3 = {—a, 1} or {l,a} with a > 0. Using 
similar arguments as in Case 1 we may assume that B\ = {—a, 0, 1} and B2 = {—a + 1, . . . , —1} in 
the first case and B\ = {0, 1, a}, B2 = {2, . . . , a — 2} in the second case (the error being sectorial). 
This gives the following possibilities for A[B\, B2,B 3 , 1] in the dimension reduced form (acting on 
M 3 ): 

Tif(yi,y 2 ,y 3 ) = 1(2/1 = a - 1,2/2 = 0)(/(a - 1,2/3 + 1) - /(a, 2/3)) 

+ l(y 2 = a - l,y 3 = 0)(f( yi ,a - 1 

r 2 /(yi,y 2 ,2/3) = 1(2/1 = 0,2/2 = a - l)(/(a - 1,2/3) - f{a,y 3 )) 

+ 1(2/2 = 0, 2/3 = a - l)(/(2/a - 1, a - 1) - /(2/1, a)) 



f{yi,a)) 



(155) 
(156) 



Since |a| < by the definition of we may further simplify these operators on Ai^: 



?l/(2/i,2/2,2/3) 
72/(2/1,2/2,2/3) 



1(2/1 = a - 1, 2/2 = 0)(/(0, 2/3 + 1) - /(0, 2/3)) 
1(2/2 = 0, 2/3 = a - l)(/(2/l - 1, 0) - /(2/1, 0)) 



(157) 
(158) 



Finally, we note that by (43) we have that \xf ' {p) is equal to the sum of the second order 
coefficients in the orthonormal representation of the microscopic flux function W. By (21) we can 
write W as 

W = s &^(y)r(y, T- Z rj)(x + X 1/2 Vz+y(l -p)~ PX 1/2 Vz - XV{z,z+y}) (159) 

y zeZy 

where \Z y \ = \y\. The sum of the coefficients of the second order terms come from the constant, 
first order and second order coefficients in the representation of the rate functions r(y,rj) and it 

gives -xJ2 y y c y,® + X 1/2 {l-'2p)J2y,\A\=iy c y,A + xY,y,\A\=2yCy,A- Recalling the definition of K from 
(108) we get that this is exactly equal to — x c + which shows that c+ = —j"{p)/2. 
Case 3: |Bi| — | B3 1 = 0. We will need to consider two cases: = 1 or 2. We first deal with 
the |-Bi| = 1 terms. Note that in that case I-B2I > 1. Using the same arguments as in the previous 
case we can show that A[B\, B 2 , {y}, y] is equal to the sum of 2/^4[{0}, {1}, {1}, 1], some terms with 
\Bi 



I-B3I = 2 and sectorial terms. Note that ^4[{0}, {1}, {1}, 1] is sectorial by Lemma 14 



only have to describe the terms with |£>i| = |£>3| = 2. We only need to do this on Ai 



dimension reduced forms of these terms vanish on M2 by Lemma 

If I Si I = I-B2I = 2 then by the previous arguments we may rewrite ^4[-Bi, B 2 ,B 3 , y] as the linear 
combination of terms A[i?i, B 2 ,B 3 , 1] with: 

Bi = B 3 = {l,a}, B 2 = [0,a]\B 1 , a > 0, or (160) 



16 



ay 



so we 
as the 



B 3 = {-a,l}, B 2 = [-a,l]\B 1 , a > 0. 



161) 



and some sectorial terms. The dimension reduced form of the operators acting 011^3 are easy to 
write down: 

Ii/(2/i,2/2) = 1(2/1 = a - l)(f(a,y 2 ) - f{a - 1,2/2)) + 1(2/1 > 0,2/2 = a - l)(/(yi - I, a) - f(yi,a- 1 

(162) 

T 2 f( yi ,y 2 ) = 1(2/1 = o)(/(o -1,2/2 + 1)- /(a, 2/2)) + 1(2/2 = a)(f( Vl , a - 1) - /(2/1, a)) (163) 
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Moreover, these can be further simplified if we restrict ourselves to ■ 

T 1 f(y 1 ,y 2 ) = > 0,y 2 =a- l)(/(yi - 1,0) - f(y u 0)) 
T 2 f( yi ,y 2 ) = = o)(/(0,itt + 1) - /(0, y 2 )) 



(164) 
(165) 



Case 4: |Bi| — | B3 1 = —1. The contribution of these terms give the part of the operator A which 
lowers the degree by one: A4 n —> At n -l- Since A* = —A from this it follows that this will be 
minus the adjoint of the part which raises the degree by one (which was discussed in Case 2). From 
this it is not hard to check that we can rewrite the contribution of these terms as c+A_ (with c+ 
from (154), and j4_ being the A- part of the TASEP), building blocks with \B\ 

2,|5 3 | 



sectorial terms. However, by Lemma 16 the terms with \B\\ 



2, 1-B3I =3 and 
3 will vanish after the 



dimension reduction so we are only left with c+^4_. 

Case 5: |Bi| — IB3 1 = —2. Since \B\\ < 3 and I-B2I > 2 the contribution of these terms is sectorial. 

Collecting all the cases we can see that Lemma [5] is proved, we just have to show that the 
contribution of the \Bi\ = 3, \B 2 \ = 2 and \Bx\ = \B 3 \ = 2 terms (see ( p7| >, ( p8| ), ( fl64| ), and gg5| |) 
is a linear combination of the operators (73) and (72) plus some sectorial terms. 



We first look at the contribution of the I B 



B, 



by proving that if we define T\ and T 2 as (164) and (165) with a 
sectorial. From the variational formula 



2 terms from (164) and (165). We will start 

T\ and T 2 — T 2 are 



1 then T\ 



(((Ti-Ti)/,^)- 1 ^!-^)/)) 



sup 2{{g,(Tx 
g£M 3 



sup 2 {a{yi,a - 1 

9&M3 yi> Q 



yi >o 



ro/))-^^-^)) 

5(2/i,0))(/(yi-l,0)-/( 2/1 ,0)) 
l,0)-/( yi ,0)) 2 <2C((/,(-5 )/)). 



((g,(-S )g)) (166) 
(167) 



The first inequality follows from ((g, (—So)g)) > C J2yi>o(9(uii a ~ 1) — 5(7/1, 0)) 2 which is a con- 
sequence of the Cauchy-Schwarz inequality (note that C depends only on \a\ < K). The same 
argument works for T 2 — T 2 . We will further modify T\ and show that it can be well approximated 
by 



f[f( yi ,y 2 ) = l(y 2 = 0) (f( yi , 0) - f( Vl 
(i.e. the difference is sectorial). This follows a similar argument: 



1,0)) 



(168) 



(((fx-ftif^-SoTHfx-fOf)) 



sup 2({g, (Ti 



sup 2 ^ «/(ltt,0)(/(ift-l,0) 



sup 2^2(g(y 1 ,0)-g(y 1 + l,0)){f(y l + l,0) 
g£M 3 yi>0 



f[)f)) - ((g, (S )g)) 
2f( yi ,0) + f( yi + l,0)) - ((g,(-S )g)) 
f(y 1 ,0))-((g,(-S )g)) 



<C£(/( yi + l,0) 
yi>o 



/( 2/1 ,0)) 2 <2C((/,(-5 )/)). 



(169) 
(170) 
(171) 



We used /(0,0) = /(1,0) = /(2,0) which follows from / G . Now note that f[ - T 2 is 

exactly the operator (72) without the constant C. Thus we just have to prove that the sum of the 
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coefficients of the operators T\ from ( 164 ) is minus the sum of the coefficients of the operators T 2 
from (165). 



This requires a bit of a book-keeping, using the arguments preceding Lemma 10 one can show 
that we only get \B\\ = \B%\ = 2 terms when |A| = 1, 2 or 3 and we have the following contributions: 
Case 1. A = {a}, a / y, 0: The contribution is y /x(n 2 —2)c yi \A[{qi, qi+1}, {<?i+2, . . . ,52}, {y, a},y] 
where qi = min(0, y, a), q2 = max(0, y, a) 
Case 2. A = {01,02}, a, 7^ y,0: The contribution is 

Ef=i 2 \/XKCy,AA[{q 1)i , {gi,i+2, . . • , Q2,i}, {y, <h}, y] with q lti = min(0, y, Oj), g 2 i = max(0, y, a* 

Case 3. A = {01,02,03}, Oj 7^ y,0: The contribution is 

Ef=i Vxc y ,AA[{q lt i, qx,i+l}, {qi,i+2, q 2;i }, {y, aj, y] with g M = min(0, y, a*), q 2yi = max(0, y, a;) 

One can also check by the preceding computations that A[B\,B2, {y, a}, y] is equal to aT[ — f3T 2 
plus sectorial terms and a — (5 = \a\ — \a — y\. We will show that in each of the three cases listed 
above the sum of the a coefficients is same as the sum of the j3 coefficients this will prove aT[ — f3T 2 
is a constant multiple of the operator (72). In order to do this we just have to prove that for any 
fixed k we have 



J2 vEf 

y,\A\=k aeA 



2/1) = 0. 



(172) 



But miraculously, this follows easily from the divergence condition (14), similarly to the last argu- 
ment in the proof of Lemma 19). This completes the first part of the proof. 

What left is to show that the contribution of the \B\\ = 3, 1-B2I = 2 terms (see (157), (158)) is a 
constant multiple of the operator (73) plus some sectorial terms. This can be done similarly as the 
\Bi 



Bo 



2 case, we will leave the details to the reader. 



□ 



Proof of Lemma [6| The proof is very similar to the proof of Lemma [5} We start with the represen- 
tation A = b 2 B 3 ,yA[Bx, B 2 , -B3, y] and note that since we are in two dimensions we only need 
to consider the terms with |JBi| < 2, the others being sectorial by Lemma 12 We will only need 
to worry about the parts of the generator which leave the degree the same or raise it, the other 
parts can be computed from the skew adjoint property A* = —A. This means that we only have 



\B 3 \ =2. 



to consider the following three cases: = 2, 1-B3I = 1, |J3i| = \B^\ = 1 and \B\ 

We start with the contribution of the terms \B\\ = 2, 1-B3I = 1. Using the same argument as in 
Case 2 of the proof of Lemma[5]we may assume that B\ = {0, y}, B 2 = and B% = y (the difference 
being sectorial). Then x 1 ^ 2 A[B\, B 2 , B3, y] is the A + part of an exclusion with jump y (and rate 
1). Using the fact the A + of a mean zero exclusion process is sectorial (see the arguments following 
( |153 )) we may replace this with the A + from an exclusion with nearest neighbor jump rates with 
drift y (again, the difference is sectorial). The total contribution of these operators will also going 
to be the A + of a nearest neighbor exclusion process and the drift is given by (154) (where y now 
runs through Z 2 vectors). 

The |-E?i| = I S3 1 = 1 terms can be treated as the similar terms in the one dimensional case: 
A[Bi, B 2 , B3, y] can be replaced with linear combinations of ^4[{0}, {e^}, {e^}, ej] (i = 1,2), terms 
with I Si I = I S3 1 = 2 and sectorial terms. The operator ^4[{0}, {e^}, {ej}, e*] is sectoria l by Lemma 
14 Finally, the building blocks with |Si| = IS3I = 2 are sectorial on M 2 by Lemma |l6 

Thus A can be written as a linear combination of A + + A^ coming from a nearest neighbor 
asymmetric TASEP and sectorial terms which proves the lemma. □ 
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5. Proof of the hard core removal lemma 
The proofs for the one dimensional case will depend on the following lemma. 
Lemma 17 (Foldout lemma). For any G\ : I? + ->• R, G 2 : Z 4 -»• R satisfying 

Gi(xi,x 2 ) = gi(xi - x 2 )l(xi = or x 2 = 0), (173) 
G2(xi,X2,X3,Xi) = Q2{x\ — X 2 , X 3 )\{x\ = X4 = or X 2 = X4 = 0) (174) 
we have the bounds 

((d, (A - So) -1 ^)) < C((gt, (A - A)- 1 ^)), ((G 2 , (A - S )- 1 G 2 )) < C((g 2 , (A - A)" 1 ^)) (175) 
where g\ : Z 2 — > R, g 2 : Z 4 — > R are defined as 

gi(xi,x 2 ) = g 1 (x 1 )l(x 2 = 0) (176) 
g 2 (xi,x 2 ,x 3 ,Xi) = g 2 {xi, \x 2 \ - l(x 2 < 0))l(x 3 = x 4 = 0). (177) 
and A denotes the usual lattice Laplacian in the appropriate dimension. 

Proof. Consider the map T ■ (x,y) — > {x — y,min(x,y)) which gives a one-to-one correspondence 
between and Z x Z + . This is like folding out Z^_ into Z x Z+. It is Lipschitz, in particular, 
nearest neighbor points are mapped to pairs with l\ distance of 1 or 2. Let q : Z+ — > R be an 
arbitrary (compactly supported) function. Define q : Z 2 — > R as 

-/ \ jq(T-Hx,y)) ify>0, 

q(x,y) = i (178) 
[q(x,-l-y) if y < 0. 

Then it is easy to check that ((Gi,q)} = ^((gi,q)}, ((q, (A — So)q)) > C((q, (A — A)q}}, and the first 
part of (175) follows from the variational formula. The second part of (175) is similar. We start 
by the 4-dimensional version of the 'folding-out' trick. Let q : Z 4 ^ — > R be an arbitrary (local) 
function. Define q : Z 4 — > R by 

( <7(T -1 (xi,x 2 ),x 3 ,x 4 ) if x 2 > 0,x 3 > 0,x A > 
q{xi,x 2 ,x 3 ,X4) = < (179) 
[ q(x\, \x 2 \ — 1(x 2 < 0), |a?3| — 1(23 < 0), \x^\ — l(x 4 < 0)) otherwise 

Then ((G 2 ,q)} = ^((g2,q)}, ((q, (A — Si)q}} > C((q, (A — A)q)) and the proof again follows from the 
variational formula. □ 

Proof of Lemma [?| In order to prove ( |76[ ) note that from the definition of F we get 

(F, (A - A)F) = 4A £ /(x) 2 + 2 £ (/(*) - f(x + e,)) 2 (180) 

where ei, e 2 are the usual unit vectors. For (/, (A — So)f) we get 

A E + 2 £ ^ " ^ x + e ^ 2 +2 ^ (/(*) - /(* + e i - e 2)) 2 (181) 

The last term can be bounded by the second one using Cauchy-Schwarz which shows that (/, (A — 
S )f) < C(F, (A - A)F). Now (F, (A - A)F) can be bounded by the right hand side of @ by 
simple Fourier computations. 



DIFFUSIVITY OF LATTICE GASES 



34 



To prove (77) we consider the following modification of Atasep- 

Af(y u y 2 ) = l( yi = 0, y 2 > 0)(/(y 2 ) - f(y 2 - 1)) + l(y 2 = 0)(/(yi) - f{yi + 1)). (182) 
The difference T = A — Atasep is given by 

Tf{y u y 2 ) = l( yi = 0,y 2 > 0)(2/(y 2 ) - f(y 2 - 1) - /(y 2 + 1)) (183) 
assuming /(l) — /(0) = 1, which we may. Using the variational formula 

(Tf, (A - 5 )- x r/) < sup 2 ( 2 f(y + 1) - /(v) - /(i/ + 2 )Mz, y) - - g(v + 1)) 2 (184) 

9 y y 

Summing by parts, it is easy to see that this is bounded by (/, (— So)f). To bound (Af,(\ — 
So)~ 1 Af) we use Lemma 17 with the function gi(x) = f(x — 1) — f{x) if x > and gi{x) = 
f(—x + 1) — f(—x) if x < 0, which gives 

(If, (A - So^Af) < C"su P 2((F(x) - F(x + l))l(y = 0),g)- {g, (A - A)g). (185) 



The right hand side is easily computed using Fourier transform as 



c 



\F(h)\ 2 \^ 



II 



_ n J_ n A + le** 1 - II 2 + |e ii2 - II 



r dt!dt 2 < C" / |F(t)| 



Of 



(186) 



which proves (77). 



Finally we prove (78) where A is defined in (71). We will assume C = 1. It is easy to see that 
Af( yi ,y 2 ,y 3 ) = l(y 2 = 0)A+f( yi ,y 3 ) + lL(y 2 = y 3 = 0)(2/(yi + 1) - /(yi) - f( Vl + 2)). (187) 



By the argument used in (184), the H—\ norm of the second part is bounded by (/, (—5b)/). From 
the variational formula, 

(1(2/2 = 0)A + /(y!, y 3 ), (A - 5 )- 1 ll(y 2 = 0)A + f( yi ,y 3 )) < (A+f, (A - So^A+f) (188) 

Using (76) and (77) the bound (78) now follows. □ 

Proof of Lemma [S[ The proof is similar to that of Lemma [7j but we need to handle a few more 
terms. We have (/, (X-So)f) = 4(F, (A- A)F) and (F, (X—A)F) can be bounded by the right hand 
side of (80 ) by Fourier computations. In order to get (81 ) we will need to bound (Aif, (\—So)~ 1 Aif) 
for i = 1,2,3 where A{ are defined in (72), (73) and (74). We will start with i = 3. 



We first separate A 3 f into two parts, we will bound the H_i norms separately. 

^3,i/(yi,y2,y3,y4) = Hyi = y2 = o)(/(y 3 + 1,2/4) - Siv^yi)) 

+ l(y 2 = y 3 = 0)(/(yi + 1, y 4 ) - f(yi + 2, y 4 )), 
A 3 , 2 f{yi,y 2 ,y 3 ,y4) = l(y 2 = y 3 = 0)(/(yi + l,y 4 + 1) - /(yi + l,y 4 )) 

+ l(y 3 = y 4 = 0)(/(yi,y 2 + l) -/(yi,y 2 + 2)) 



(189) 
(190) 
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We modify these a little but by introducing 



A 3 ,if(yi,y2,y3,y4) = Hvi = V2 = o)(/(y 3 + i,y 4 ) - f(y3,yd) (191) 

+ Hv2 =y 3 = 0, yi > l)(/( yi - 2, y 4 ) - f(yi - l, y 4 )) 
A 3 , 2 f(yi,y2,y3,yA) = Mm = y 3 = o)(/(yi,y 4 + 1) - /(yi,y 4 )) (192) 

+ l(y 3 = y 4 = 0,y 2 > l)(/(yi,y 2 - 2) - /(yi,y 2 - 1)) 

The fact, that ((i 3 , 4 - I 3) i)/, (A - 5c)" 1 (A 3 ,i - A 3)i )f) for i = 1, 2 are bounded by C((/, (-S )/» 
follows the same way as by the argument used in ( |184 ). 

To bound the norm of A 3 ^f we use Lemma (17) with the following setup: 



G 2 (xi, x 2 , x 3 , x 4 ) = A 2j if(x 2 ,x±, xi,x 3 ) 

f{x + 1, y) - fix, y) = Fix + 1, y) - Fix, y) 
fi-x - 2, y) - fi-x - 1, y) = F(x + 1, y) - y) 



92ix,y) 



if x > 
if x < 



(193) 



where we used / S .M 3 for the x = case. Then y 2 defined in (177) is given by 
g 2 {xi,x 2 , x 3 , x 4 ) = (F(xi + 1, x 2 ) - Fixi,x 2 ))l{ X3=X4=Q y. 



(194) 



Now (y 2 , (A — A) g 2 ) can be computed explicitly using Fourier-transform of F which leads to 



(52, (A -A)- 1 ^) <C 



f2 



t,tt]4 a + tf + q + tg + 1\ 



|F(ti,t 2 )| 2 dti(it 2 dt 3 dt 4 . 



(195) 



Integrating out t 3 and t 4 gives the first term on the right hand side of (81 ). We can bound the H-\ 
norm of A 3)2 f exactly the same way (using /(x,y) = /(y,x)). 

To bound (A 1 f, (A-S'o)" 1 ^!/) note that with /(y) = /(0, y) we have Aif{yi,y 2 ) = A TASEP fiy 1 ,y 2 ) 
and we could use the bounds from the proof of Lemma 7 The role of Fit) is played by f* n F(ti , t 2 )dt\ 



which gives the second term on the right hand side of (81 ). 

Finally, one can show that (A 2 f, i\ — So)~ 1 A 2 f) is bounded by C(A\f, i\ — So)~ 1 Aif) using the 
same argument as in the end of the proof of Lemma [7| This finishes the proof of Lemma [8| □ 

Proof of Lemma^ Consider the symmetric extension /* of / on (Z 2 ) 2 — > M. This is given by 
/*(xi,yi,x 2 ,y 2 ) = /({(xi,yi), (x 2 ,y 2 )}) if (xi,yi) / (x 2 ,y 2 ) and zero otherwise. By Lemma 4.3 
of [19] we can bound (/, (A — So)/) with C(/*, (A — A)/*) and the same estimate holds for the 
quadratic forms with ((•,•)). A simple Fourier computation shows that 



((/*,(A-A)/4 <C 



(A + \u\ )f*iu, —u)du. 



(196) 



Here /*(n, v) is the Fourier transform of /*, and it is clear that /*("U, —it) = -F(it) where F is defined 
in (g. This proves Q. 

The second part of the lemma follows directly from Lemma 3.2 and 4.5 of |19j . Note that 
although in the paper they only deal with the (a, b) = (1,0) drift, the arguments can be extended 
to the general case without any difficulty. □ 
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6. Appendix: Green-Kubo formula 



Lemma 18. Let Ca = J2yVi Er{y,rf). Then 



ft* e s pt f s 

tD ii (t) = C i , i t + 2 X ((w^^^Wiis^duds- {{v i {s),e uL *v l {s)))duds, (197) 

Jo Jo Jo Jo 



Proof of Lemma 18 We will only show the d = 1 case, the general case is similar. The first part 



of the proof is standard. We start with the martingales 

M x (t) = r] x (t) - n x (0) - [ ^2 ( r (y> Tx-yVhx-yil - Vx) - r(y, T x rj)i] x (l - rj x+y )) ds (198) 



o 



y 

Using the notation V z f(x) = f(x + z) — f(x) and the identity J2x(^ ' zo)b = ^ x (V_ 2 6)a we get 
*£x 2 ( Vx (s), Vo (0)-p) = J2 x2 (Vx(0) + M x (t), m (0)-p) 



X X 

rt 



+ ^2 x2 (f ^2(^-y( r (y, T xV)Vx(l -rj x+y )))ds,rio(0) - p). 
x ® y 

The first term vanishes by the independence of {i]x(0)} under tt p and the martingale property. Now 
we use the generator L* of the time reversed process to get 

r]x (0)= Vx (s)+ f £ (r(^r a7 7^ + ,(l-r/:)-r( 2j r^>*(l-r7*_J)d S + M*( S ) (199) 
\z\<K 

where r]*{u) is the time reversed lattice gas process started from r/*(0) = rj(s). This gives 

5> 2 Ms), 770(0) - p) = f ((-2xy + y 2 ) (r(y, 77)7/0(1 - Vy )) ,r fx -p + M*(s))ds (200) 
~ t~r. Jo 



x x,y 



+ E / / ( (~ 2x y + y 2 ) ( r (y> v)vo0- - %)) , - E v -* ( r ( z > - €)) ) duds 

x,V J ° J ° \ \z\<K I 

(201) 



In the last term n x = r] x (s) and 77* = n*(u) and V_ 2 acts in the x variable. We will compute ( 200[ ) 



and (201) separately. 



We start with (201). By moving the gradient V_ 2 onto the first term we get 



^2yz [ f (r(y,rj)rio(l-%)-,r(z,T x ri*)r)* +z (l-ril))duds 

x,y,z J J 

= 2 f [\(£yr(y,v)rio{l-Vy),Y, zr ( z > r >*K(l-V*a))duds (202) 

JO J0 y 

This gives 

X - 1 J2 x2 (Vx(s),Vo(0)- P ) = Ct + 2 X - 1 f f S ((W(s),e uL *W*(s)))duds (203) 

Using (a,b) = ((o+6)/2, (o+6)/2)-((o-6)/2, (o-6)/2) and introducing v := W-W*,W := W+W* 
we can rewrite the integral term as 

ij p s pi ps pi ps 

/ ((W{s),e uL *W*(s)))duds = / / ((W(s), e uL *W(s)))duds- / / ((v(s), e uL * v(s)))duds 
Jo Jo Jo Jo Jo 
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Recall that w = W - V{W - (W). By the conservation of £77 we have ((ViW, e uL * ViW (s))) 
((VxW,ViW(s)}) and ((u>, e uL *7W(s))) = 0. This gives 



t rs 



t rs 



{{W(s),e uL W{s)))duds = 2 / {(w(s),e u ^ w{s)))duds + t z ((V l W,T 1 W)) (204) 
JO Jo JO 

In order to compute the last term we introduce rj\ = YlxeAVx- If |A| = k then 

tja = Y\ (x 1 ^ 2 Vx + p) = p k + X l ^ 2 p k ~ l E ^ x higher order terms. (205) 

zeA rreA 

Since W is a polynomial we can write W = J2 a AVA- Then V\W = x 1 ^ 2 J2a ca/»' A ' _1 J2x.eA Vx and 



fW,^ = * 2 «Ep |am E ^.E a AP |A|_1 E *» = * 2 (E^ |Ahl i A i 

A xeA A zeA \ A / 



Note that j(p) = EW = EW = £ A a A p |A| which means that 



(206) 



(207) 



(208) 



The computation of (200) relies on a bit of a miracle, which we state as a lemma. 
Lemma 19. 

x _1 E / ((-^y + y 2 ) {r(y,v)M^ - Vy)) ,Vx - p + m;(s)) = tY / y 2 (r(y,v)) ■ 

x,y y 

Proof. By stationarity, and because M* is a martingale the sum on the right is equal to 
*X "TE^" 2 ^ + y 2 )r(y,ri)r/ (l -%)),% -p). 

2/ as 

If we consider the x = 0, x = y terms then we get 

'x -1 E y 2 ( r (y> ^0(1 - vy)(vo - vy)) = t E y 2 ( r (y, v)) 

y y 
For the other terms, by the independence of 77$,% when i / j we have 

tx~ 1/2 E E (~ 2x y + y 2 )( r (y> rivoQ- - = *x 1/2 E< r (^ v),J2(~ 2x y + y 2 )^>- ( 21 °) 

!/ x^0,y y x 

We just need to show that the term on the right is zero. Let V\ be the projection to M\ and 
denote Vir(y, 77) = Y,x Cy^Vx- (We can assume that Cy y y — c^o — 0.) Then 

^2(r(y,v),^2(-^y+y 2 )vx) = EE^fe 2 - 2 ^) = ^2 c vA(y- 

y x y x y,x 

= E a E c v> x ~ E a E c v> x 

a \y— x\=a a y,\x\=a 



(209) 



x) 2 -x 2 ) 



(211) 



By the divergence condition (14) the sum ^}2yVir{y,ri)(riy — tjq) must also be a gradient. By 
inspecting the coefficients of fjy, y + a we get that this is equivalent to 



E 

\y-x\=a 



-y,x 



E C V> 2 



for all a. 



y,\x\=a 



This shows that (211) and thus (210) will vanish. 



(212) 
□ 
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